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1  i^XEIDING 


De  Universele  Dynamometerbank  (UDB)  is  een  testbank  voor  voeituigen  met  een  massa  tot 
100000  kg.  Op  de  testbank  kunnen  rijsimulaties  uitgevoerd  worden.  De  electromagnetische 
motoren  leveien  een  snelheidsafhankelijk  koppel  M  op  de  aandrijfassen  van  het  voertuig  dat 
gegeven  wordt  door  een  vieital  co^fficienten.  de  z.g.  K-factoien: 

(1)  M  =  12  K,  +  (3/50)  KjO)  +  (3/1000)  KjOp  +  (5/2)  K4(d’  (to  de  hodcsnelheid) 

Een  probleem  vonnt  de  bepaling  van  geschikte  instelwaarden  voor  de  K-factoren  Kj,  K2  en  K3, 
zodat  de  UDB  een  realistisch  koppel  leveit.  K4  kan  eenvoudig  berekend  worden  als  het  product 
van  de  voertuigmassa.  het  kwadraat  van  de  dynamische  wielradius  en  een  schaalfactor.  Kj,  K2  en 
K3  coiTCsponderen  met  wrijvingscoiiflficientea  Kj  ai  K2  conespondercn  met  de  rolweerstand. 
Deze  is  in  weikelijkheid  afhankelijk  van  factoren  zoals  de  samenstelling  van  de  rijbaan,  de 
bandenspanning,  de  temperatuur  de  bandensooit  en  de  snelheid.  Op  empirische  wijze  is  bepaald 
dat  de  rolweerstand  lineair  toeneemt  met  de  snelheid.  Kj  correspondeeit  met  de  rolweerstand 
onafhankelijk  van  de  snelheid.  K2  correspondeeit  met  de  lineair  snelheidsafhankelijke 
rolweerstand.  K3  correspondeeit  met  de  luchtweerstandscoiiflici&iL 

Voor  de  bepaling  van  de  wrijvingsco(!fficienten  is  een  bewegingsvergelijking  q}gesteld  wdke  het 
uitrollen  van  een  voertuig  modelleert.  De  oplossing  (0(0  van  deze  diffeientiaalveigelijking 
beschrijft  dan  de  uitrolcurve.  K],  K2  en  K3  heten  optimaal  bepaald  te  zijn  als  de 
differaitiaalvergelijking,  waaiin  de  waaiden  van  Kj,  K2  en  K3  gesubstitueerd  zijn.  een  cqdossing 
heeft  die  zo  goed  als  mogelijk  overeoikomt  met  de  gemeten  uitrolcuive. 

De  optimak  K],  K2  en  K3  worden  bepaald  door  een  noimfiinctie  te  optimalizeren  met  een 
optimalizeringsmetfiode. 

bi  het  documeitt  zal  eerst  het  fyrisdie  model  beachieven  wtmlen.  De  resultaten  vn  het  ondenoek 
aan  de  difTemttiaalvergelijking  wordmi  gqrreaeniead  en  een  iwim  wcHdt  gekozen  ten  opaidtte 
waarvan  geoptimalizeerd  moet  gam  worden.  Afhankdi|c  van  de  benadering  van  het  fyiiscfae 
probleem  ktmnen  veiscdiillende  nonnfuncties  opgesidd  wordea  UH  deze  normftmcties  wordt  een 


TNO-rapport 


Pagina 

7 


keuze  gemaakt.  Vervolgens  worden  twee  optimalizeringsmethoden  besptoken.  waarvan  er  66n 
gekozen  wordL 

Tot  slot  wotden  de  resultaten  van  de  optimalizeringsmethode,  toegepast  op  eea  vijftal  gemeten 
uitrolcurves  besptoken. 

Hieronder  volgen  enkele  convendes,  notaties  en  afkoitingen.  Een  vet  gedrukte  letter  geeft  een 
verzameling  aan.  Tekst  tussen  hoekige  haken  geeft  een  eenheid  aan. 

R  verzameling  van  ie£Ie  getallen 

M-dimensionale  verzameling  van  regie  getallen 
D  c  RM  D  is  een  deelverzameling  van  RM 

A  de  effectieve  oppervlakte  van  het  voeituig  [m^] 

g  de  gravitatie  constante:  9,8  m  / 

m  de  massa  van  het  voeituig  [kg] 

r  de  dynamische  wielradius  van  een  wiel  van  het  voeituig  [m] 

p  de  soortelijke  massa  van  lucht:  1 ,23  kg  /  m^ 

<0  de  hoeksnelheid  van  een  wiel  van  het  voeituig  [rad  /  s] 

(O'  de  eeiste  afgeleide  van  de  hoeksnelheid  naar  de  tijd  [rad  /  s^] 
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2  HET  FYSISCHE  MODEL 


Het  fysische  model  is  de  bewegingsvergelijking  van  het  voertuig.  Uit  de  bewegingsvergelijking 
en  een  verzameling  meetwaaiden  Icunnen  de  K-factoren  bepaald  woiden. 

Er  wordt  verondersteld  dat  een  verzameling  W™  met  N  meetpunten  beschikbaar  is: 

(2)  W">={(ti.a)«»i)IO<ti.,<ti<tMax,i  =  2,3,....N-l;to  =  0; 

0  S  (D"i  ^  ay"Max,  i  =  0. 1, ....  N  - 1;  oa^o  >  0 ) 

waaibij  (tj,  oy"i)  een  meeqnint  is  en  waarbij: 

tj  het  tijdstip  van  de  i*'^  meting  in  seconden  na  de  eerste  meting. 

(lyOj  de  gemeten  hoeksnelheid  van  een  wiel  in  radialen  per  seconde  op  tijdstip  tj. 

Opgemerkt  zij  dat: 

zonder  beperking  der  algemeenheid  aangenomen  wordt  dat  de  eerste  meting  gemeten  is 
op  tijdstip  0  seconden;  dus  to  =  0. 

om  te  kunnen  spieken  van  een  uitrolcurve,  (0’"o  positief  moet  zijn. 


De  bewegingsvergelijking  is  een  eerste-orde  difTerentiaalvergelijking  zonder  randvoorwaaide.  De 
bewegingsvergelijking  kan  op  verschillende  wijzen  geschieven  wotden. 
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De  fysische  bewegingsvergelijking  met  f],  f2  en  c^: 

(3*)  m  r2  ©'  +  m  g  r  fj  +  m  g  r2  fj  m  +  A  p  r3  0)2  =  0 
waaibij; 


f]  de  snelheidsonafhankelijke  wrijvingscogfficiSnt  voorstelt 

f2  de  lineair  snelheidsafhankelijke  wrijvingscoefiicitint  voorstelt 

Cy,  de  kwadratisch  snelheidsafhankelijke  luchtweerstandscodfficignt  voorstelt 

Meik  op  dat  helling,  de  windsnelheid  en  -richting  niet  meegenomen  woiden  in  het  fysische 
model. 


De  mathematische  bev^egingsvergelijking  met  a,  b  en  c: 
(4)  0)'sa  0)2-»-b  0)-t-c 


De  coeflicienten  a.  b  en  c  woiden  mathematische  co(!fficienten  genoemd.  Vanuit  mathematisch 
standpunt  bezien  is  de  schiijfwijze  (4)  te  prefereren  aangezien  deze  het  meest  compact  is. 

Voor  de  inteipretatie  van  de  resultaten  kunnen  de  mathematische  co^iflicientcn  het  beste 
omgerekend  woiden  naar  fysische  coefficiSnten.  Voor  het  maken  van  de  instellingen  aan  de  UDB 
kunnen  a.  b  en  c  omgerekend  wordra  naar  K|,  K2  en  K3. 


Val|MH  KnATriiilf^  ii  an  «un  Hoar  1/2  «w  de  fcMkelyk.  Dit  bedt  M  (cvolg  da  aBe  (^•waadai  im 

ditdocwnanimal  wadgaddiodmoeiaiwoideaonediaiecbeweiidaiayeifciiitea 
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De  formules  om  K-factoren  en  fysische  co£fficienten  te  berekenen  uit  de  mathematische 
coSfiiciiinten  volgen  hieionder: 


m 

r 

- c 

12 

f,  =- 

—  c 

g 

50  m 

u 

1 

3 

h  -  * 

g 

1000  m 

m 

3  Apr 


Een  oplossing  van  een  differentiaalvergelijking  is  niet  uniek  bepaald  zonder  randvoorwaarde. 
Derhalve  wordt  verlangd  dat  de  oplossing  op  tijdstip  to  een  (onbekende)  hoeksnelheid  (Dq 
aanneemt; 

f  (O'  =a  co^  +  b  co  +  c 
(5)  \ 

[  aKto)  =(0t, 


Meric  op  dat  bet  aantal  onbekende  codfficidnten  nu  met  6€n  is  toegenomen  ! 

De  'beste'  K-factoren  corresponderen  met  die  mathematische  codfficidnten  waaivcwr  de  oplossing 
van  stelsel  (S)  het  'beste'  past  bij  de  meetwaaidenverzameling  Wi". 
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3  ONDERZOEK  AAN  DE  DIFFERENTIAALVERGELUKING 


Het  eeiste  doel  is  het  opstellen  van  een  nonnfunctie,  die  geoptimalizeeid  de  'beste' 
wrijvingscodfficidnten  oplevert.  Daarvoor  is  onderzoek  van  differaitiaalvergelijking  (5) 
noodzakelijk. 

Het  blijkt  mogelijk  te  zijn  de  differentiaalvergelijking  (5)  op  te  lessen.  Van  gnx)t  belang  voor  het 
toepassen  van  optimalizeringsmethoden  zijn: 
existentie  van  oplossingen 
uniciteit  van  een  oplossing 
continuVteit  van  een  oplossing 

differentieerbaarheid  van  een  oplossing  voor  snelle  conveigentie 


3.1  Existentie  en  uniciteit  van  een  oplossing 


M.b.v.  een  tweetal  stellingen  is  te  bewijzen  dat  de  differentiaalvergelijking: 

[  U’(t)  =  f(t,  U(t)) 

(6)  \ 

lu(a)  =Co 


een  oplossing  heeft  en  dat  deze  oplossing  uniek  is. 


Aanname  D  =  {(t,  x)  I  a  S  t  S  a  +  6,  llx  -  cqII  S  y},  en  f :  D  RM  is  continu  (5.  y>  0  en  ll.ll 
een  willekeurige  norm). 

SteUing  Onder  aaimame  I  is  er  een  getal  ^  met  a  <  p  ^  a  •«-  8,  zo  dat  (6)  een  oplossing  U 
heeft  op  het  interval  [a,  p]. 
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Aanname  D  c  en  f :  D  -»  RM,  zd  dat  bij  elke  begrensde  B  c  D  een  L  (=  Lg)  bestaat 
met  llf(t,  x)  -  f(t ,  2011  ^  L  llx  -  jll  (voor  alle  (t,  x),  (t,  2)  e  B). 

Stelling  Onder  aanname  II  bestaat  er  voor  (6)  hoogstens  ddn  opiossing. 


Met  f(t,  x)  =  a  x^  +  b  X  +  c  is  bet  evident  dat  aan  aanname  I  voldaan  kan  worden.  Er  bestaat  dus 
volgens  stelling  I  tenminste  ddn  opiossing  voor  (6). 

Met  f(t,  x)  =  a  x^  -«■  b  X  +  c  is  00k  voldaan  aan  aanname  II;  omdat  B  begrensd  is,  is  ook  L 
begrensd. 

Er  bestaat  dus  volgens  stelling  II  hoogstens  66n  opiossing  voor  (6). 


Dus  als  a,  b,  c  en  (Oq  gegeven  zijn,  dan  bestaat  er  precies  6dn  opiossing  voor  (6). 


3  J  Continui'teit  en  diiTerentieerbaarheid  van  de  opiossingen 


StellinglOI; 

Zij  gegeven  het  stelsel  differentiaalveigelijkingen  (6)  met  f :  [0,^]  *  D  -» R^,  en 
f  e  CPX[a,p]  *  D)  (=  verzameling  van  alle  p  maal  continu  diffetendeeibare  ftmcties  op 
[ot,p]  ♦  D),  dan  is  dc  o{dossing  U(t)  van  (6)  een  element  van 

Gevolg:  als  a,  b,  c  en  o\)  gegeven  zijn,  dan  is  de  bijbehorende  opiossing  (D(t;  a,  b,  c,  (Dq)  oneindig 
vaak  continu  differentieeibaar  op  een  zeker  interval  [a,p].  Er  geldt  immers  voor 
f(t,  x)  =  a  x2  +  b  X  +  c  dat  f  e  C"l({ot.P]  *  R. 
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3J  Opiossingen 

De  diffeientiaalvergelijking  (5)  is  oplosbaar,  d.w.z.  voor  vaste  coefficienten  a,  b,  c  en  (Dq  kan  cr 
een  expliciete  fiinctie  ©(t;  a,  b,  c.  ©o)  gevonden  worden,  die  voldoct  aan  de  volgende  eisen: 

(i)  ©  is  gedefinieeid  op  een  interval  J  c  R. 

(ii)  V  t  e  J  is  (t,  ©(t;  a,  b,  c,  ©o))  e  [deiinitiegebied  f(t,  x)  =  a  +  b  x  +  c] 

(iii)  V t e  J geldt:  ©'  =  a©2  +  b©  +  c 

(iv)  ©(to;  a,  b.  c.  ©o)  =  ©0 


De  differentiaalvergelijking  (5)  is  oplosbaar  m.b.v.  separatie  der  variabelen.  De  oplossing  van  (5) 
blijkt  zich  op  te  splitsen  in  drie  (deel)oplossingen  ©j,  ©2  en  ©3.  De  opiossingen  zijn  gecontroleerd 
door  substitute  van  de  opiossingen  in  de  differentiaalvergelijking  (5)  en  door  verificatie  van 
®.(to)  =  (Oq,  zowel  met  de  hand  als  met  het  symbolische  mathematische  programma  'Mathematica'. 

De  opiossingen  van  difTerentiaalvergelijking  (5)  zijn: 


2  a  ©0  Tan [ 


2 
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-(b2  t  +  4  a  <0o  +  2  a  b  t  (Oq) 

(02[t;  a,  b,  c,  (Oq]  =  - 

2a  (-2  +  b  t  +  2  a  t  ©o) 

Sqrt [b^  -  4  a  c]  t 

©3[t;  a,  b,  c,  ©q]  =  -((-2  c  +  2  E  c  -  b  ©o  + 


Sqrt[b2  -  4 

a 

c]  t 

E 

b  ©0  +  Sqrt [b^  -  4  a  c]  ©o  + 

Sqrt (b^  -  4 

a 

c]  t 

E 

Sqrtlb^  -  4  a  c]  ©q)  / 

Sqrt [b2 

- 

4  a  c]  t 

(-b  +  E 

b  -  Sqrt[b^  -  4  a  c]  - 

Sqrt [b^  -  4 

a 

c]  t 

E  Sqrt [b^  -4a  c]  -2a  ©o+ 

Sqrt [b^  -  4  a  c]  t 
2  E  a  ©o)) 


Er  geldt: 


(7)  ©[t;  a,  b,  c,  (Oq] 


©I  [U  a,  b,  c,  (Oq]  als  -  4  a  c  <  0 
©2[t;  a,  b,  c.  ©o)  als  b^  -  4  a  c  =  0 
©3[t;  a,  b,  c,  ©oJ  alsb2-4ac>0 


3.4  Limieten 


De  afzondeilijke  opiossingen  zijn  continu.  ©  daarentegen  behoeft  niet  continu  te  zijn  als  fimctie 
van  a,  b,  c  en  Q\).  Het  oveiigangsgebied  van  de  opiossingen  ©],  en  ©3  dient  daaitoe  onderzocht 
te  woiden. 

Het  overgangsgebied  van  de  opiossingen  is  de  volgende  verzameling  piinten: 


( (a.  b,  c,  ©o)  e  I  b2  -  4  a  c  s  0 } 
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Het  gedrag  van  de  oplossingen  in  het  overgangsgebied  is  ook  numeriek  van  belang.  De 
oplossingen  kunnen  slecht  geconditioneerd  zijn  in  het  overgangsgebied,  d.w.z.  mathematisch 
gezien  welgedefinieeid,  maar  indien  geSvalueerd  met  floating-point  operaties  in  absolute  waaide 
extreem  groot. 

In  oplossing  tOi  wordt  voor  Tan(x)  de  lineaire  teim  van  de  Taylorreeks  van  Tan(x)  gesubstitueerd. 
Substitutie  is  toegestaan  omdat  in  het  overgangsgebied  x  =  (b2-4ac)*t/2-»0. 

Oplossing  (Oj  herschreven  na  substitutie  van  x  voor  Tan(x): 

(  Sqrt  [-b^  +  4ac]  (2act-2acto  +  2atDo  +  ®l^t(Oo-abto 
Wq))  / 

(  2  a  Sqrt[-b2  +  4  a  c]  -  a  b  Sqrt[-b^  +  4  a  c]  t  + 

a  b  Sqrt[-b2  +  4  a  c]  tg  -  2  a^  Sqrt[-b2  +  4  a  c]  t  ©g  + 

2  a^  Sqrtt-b^  +  4  a  c]  tg  Wg) 

De  bovenstaande  formule  met  b^  /  (4  a)  gesubstitueerd  voor  c  blijkt  identiek  te  zijn  aan  oplossing 
(02. 

Op  soortgelijke  wijze  kan  aangetoond  woiden  dat  (03  identiek  is  aan  oplossing  (O2,  door  1  +  x  te 
substitueren  voor  e*  (eerste-orde  Taylorreeks). 

Er  geldt  dus: 

(8)  lim  (D(t;  a,  b.  c,  OQo)  =  (a2(t;  a,  b.  c.  ©0) 

Y 

waarbij  y  een  willekeurig  pad  is  naar  een  punt  waarvoor  b^  -  4  a  c « 0  geldt  (8)  geldt  voor 
willekeurige  a,  b,  c.  ©q.  De  (x>ntinuneit  in  het  overgangsgebied  b^  -  4  a  c  =  0  is  dus  gegaranteeid. 
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33  Particle  afgeleiden 


Optimalizeringsmethoden  maken  veelal  gebruik  van  paniOe  afgeleiden  tb.v. 
convergentiesnelheid.  De  paitiSle  afgeleide  van  co  naar  a  blijkt  editer  niet  uniek  gedefinieeid  te 
zijn.  Dit  blijkt  uit  de  beiekening  van  de  pardSle  afgeleiden  van  0)101(03  naar  a  in  bet  punt  (t  =  S, 
a  =  -l,b  =  -2,  c  =  -l,o)o  =  0)(In  feite  is  de  limiet  naar  dit  punt  genomen  aangezien  o>i  en  0^  niet 
gedefinieeid  zijn  alsb2-4ac  =  0en  dus  zeker  00k  de  partible  afgeleiden  niet  gedefinieeid  zijn). 
Het  niet  uniek  gedefinieeid  zijn  van  de  parUdle  afgeleide  heeft  tot  gevolg  dat  snel  conveigerende 
optimalizeringstechnieken  die  gebruik  maken  van  partidle  afgeleiden,  uitgesloten  zijn,  aangezien 
deze  technieken  verlangen  dat  de  eerste-orde  paitiSle  afgeleiden  tenminste  (x>ntinu  zijn. 

Veider  onderzoek  aan  de  paitidle  afgeleiden  ten  behoeve  van  een  snel  conveigerende 
optimalizeringsmethode  is  derhalve  niet  zinvol. 
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4  NORM 


Om  twee  veizamelingen  co<!fnci£nten  met  elkaar  te  kunnen  veigelijken  op  'kwaliteit'  (wdke 
verzameling  co^fficii^ten  levert  een  oplossing  die  het  meeste  lijkt  op  de  meetwaardencuive)  is  de 
introductie  van  een  norm  noodzakelijk. 


Stel  o>  is  een  oplossing  behorende  bij  een  parameterverzameling  P.  Laat  O  de  verzameling  zijn 
van  alle  functies  f :  R  met  M  het  aantal  onbekende  co6fiici€nten  in  P.  Dus  o>  e  O. 


Een  voor  de  hand  liggende  keuze  voor  de  norm  is  de  som  van  de  kwadratische  verschillen: 
ll.llj  :  fl  -4  Rjo 

N-l 

©  -4  Z  (©"i  -  ©(tj.  P))2 

Andere  mogelijkheden  zijn: 

Oe  som  van  de  absolute  verschillen; 

"•H«m  :  a-* 

N-l 

©  -»  1  -  ©(tj,  P)l 

taO 

De  maximumnoim; 

H-Knm  £2 

©  -4  max  -  ©(tj,  P)l 
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Een  minder  gebruikelijke  imrm  is  de  volgende: 

IIJIp  :  Rjo  (P>0) 

N-l 

fi)  -»  X  lay“i  -  <0(1;,  P)l*^ 


Formeel  gezien  behooit  tnj  II.II2  nog  de  wortel  genomen  te  worden  en  bij  ll.llp  de  (l-fp)<>B-machts 
wortel.  Aangezien  de  norm  slechts  aangewend  woidt  voor  de  vergelijking  van  opiossingen,  spelen 
de  woitels  geen  rol. 

De  ll.llp  valt  af  als  norm  omdat  deze  noim  de  spreiding  van  de  afwijkingen  van  de  meetwaaiden 
van  de  oplossing  van  de  diffeientiaalvergelijking  niet  minimalizeert 

Omdat  veel  optimalizeringstechnieken  maken  gebmik  van  paitiSle  afgeleiden,  of  bij  bet  ontbreken 
daarvan,  van  'finite  differences'  voor  de  benaderingen  van  partidle  afgeleiden  is  bet  niet  oontinu 
diffeientieerbaar  zijn  van  H.llmn  oipraktiscb.  Er  geldt  voor  optimalizeringstecbnieken  die  gebmik 
maken  van  finite  differences  dat  des  te  meer  discontinuiteiten  in  de  panicle  afgeleidoi,  des  te 
kleiner  de  kans  dat  er  een  globaal  optimum  gevonden  wordt.  Er  is  een  goed  altematief  waardoor 
de  nadelen  van  il.Hmu  ontweken  kunnen  wordrai. 

II.II2  en  ll.llp  zotgen  voor  een  zo  klein  mogelijke  spreiding  van  de  afwijkingen  en  leiden  tot 
betzelfde  optimumlf^l.  Il.llj  verdient  de  vooikeur  boven  ll.llp  aangezien  ll.llj  minder  rekenweik 
verge 

De  keus  voor  de  norm  valt  op:  II.II2. 
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5  NORMFUNCTIE 

De  nonnfunctie  is  gedefinieerd  als: 

N-l 

F(P)  =2 

pO 

De  bepaling  van  (D(tj;  P)  legt  F(P)  vast.  Erzijn  twee  manieren  om  (o(t;  P)  te  bq>a]ai. 


5.1  Numerieke  normfunctie 

Bij  de  numerieke  nonnfunctie  woidt  geen  gebniik  gemaakt  van  de  exacte  oplossingen  van  (5), 
maar  wordt  de  functie  m  numeriek  benadeid.  Uitgangspunt  voor  de  numerieke  methode  is  het 
volgende  stelsel: 

<B'  =a<»2  +  ba)  +  c 

(9) 

l<i()(to)»(OWo 

De  randvoorwaaide  van  het  stelsel  (9)  is  gelijk  aan  de  eerste  gemeten  hoeksneDieid  van  een  wid 
van  het  voenuig.  Dit  gegeven  wordt  gebniikt  voor  de  numerieke  benadering  van  de  functie  o. 

Als  nonnfunctie  wordt  de  op  II.II2  gebaseerde  functie  F  gebniikt: 

N-I 

(10)  F(a,  b,  c) «  I  («»»i  -  «Kq;  a.  b.  c))2 

De  a.  b  en  c  waarvoor  F(a.  b,  c)  minimaal  is,  zijn  de  tamale  mathematische  codflid&itea 
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De  flUKtie  oXt;  a,  b,  c)  woidt  numeriek  benaderd  m.b.v.  de  methode  van  Euler  (Zie  bijlage  A). 

(10)  kan  verfijnd  woiden  door  stelsel  (S)  als  uitgangspunt  te  nemen  en  de  volgende  normlunctie  te 
gebruiken: 

N-l 

(11)  F(a.  b,  c.  oio)  =  I  (ay"i  -  ©(tj;  a.  b,  c,  Wo))^ 

1=0 

De  optimale  mathematische  codfficienten  zijn  die  a,  b,  c  en  (Dq.  waarvoor  F(a.  b,  c,  oao)  minimaal 
is. 

De  verfijning  heeft  tot  gevolg  dat  de  invloed  van  een  foutje  in  (D’"o  vennindeid  wordt  dooidat  de 
randwaarde  kan  varieren. 

De  fiinctiewaaiden  van  (d  kunnen  op  een  willekeurig  aantal  decimalen  precisie  berekend  worden 
m.b.v.  de  methode  van  Euler  door  een  voldoende  kleine  stapgrootte  te  nemen  en  de  beweikingen 
uit  te  laten  voeien  met  voldoende  piecisie. 


5  J  Exacte  nomifunctie 


Deze  methode  gaat  uit  van  stelsel  (S).  De  normfiinctie  is  (1 1),  waarbij  O)  bepaald  wordt  m.b.v.  de 
exacte  opiossingen  van  (5).  De  a,  b,  c  en  cog  waarvoor  F(a,  b,  c,  cOq)  minimaal  is.  zijn  de  opdmale 
a,  b,  c  en  (liV). 

Merit  op  dat  de  exacte  noimfunctie  ctmtinu  is. 


53  Kcuie  votN*  de  normftincUe 

De  exactt  normfiinctie  verdient  de  vooriteur  boven  de  numerieke  nonnfimctie  omdat  indkn 
exacte  opiosidi^en  besdAbaar  zijn  het  onzinnig  is  om  numeridte  benaderingen  uit  te  voemi.  ook 
al  leidt  het  kiezmi  voor  exacte  qplossingen  tot  meer  rdcemijd.  hetgeen  in  dit  gevai  overigens  niet 
zo  is.  Het  monmidt  benadeien  van  eo  kost  meer  rdcentijd  (0(N^)  dan  het  gebnrik  van  de  exade 
opkmingen. 
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6  OPTIMALIZERINGSMETHODEN 


Bij  de  keuze  van  de  optimalizeringsmethode  voor  het  optimalizeTen  van  de  normfunctie  is  het  van 
belang  dat  de  optimalizeringsmethode  de  opdmale  codfficidnten  moet  kunnen  bepalen  binnen 
fysisch  vastgesteide  intervallen.  De  intervallen  definidren  het  zogeheten  zoekgeUed. 

Een  tweetal  routines  woidt  beschouwd:  BCONF  en  RNLIN.  Deze  routines  zijn  afkomstig  van  de 
IMSL  library.  Beide  routines  kunnen  fiincties  optimalizeren  van  het  volgende  type: 


6.1  IMSL-routine  BCONF 

De  BCONFtC?)!  optimalizeringsmethode  maakt  gebruik  van  een  quasi-Newton  methode.  Finite 
differences  worden  gebruikt  voor  de  bepaling  van  de  paititile  afgeleiden  van  de  normfunctie. 

BCONF  biedt  de  mogelijkheid  om  expliciet  intervallen  aan  te  geven,  waarbinnen  de  codflicidnten 
moeten  liggen.  Er  woidt  niet  gegarandeerd  dat  er  binnen  het  zodtgebied  geen  andere  coSflicienten 
te  vinden  zijn  die  een  beter  optimum  hebben.  Dit  probleem  behooit  eigenlijk  ondeizocht  te 
worden  in  een  functietmdeizoek.  Vanwege  de  complexiteit  van  de  exacte  normfunctie  (m.n.  ten 
gevolge  van  de  dimensie)  en  de  data-afhankelijkheid  is  dit  analytisch  een  vrijwel  cmmogelijke 
zaak. 

De  afgewezen  numeiidte  normfunctie  biedt  als  vooideel  dat  de  normfunctie  in  feite  i^ts  anders 
is  dan  een  data-afhankelijk  4*-oide  polynoom.  Het  gevolg  hiervan  is  dat  de  numerieke 
normfunctie  geen  siiigulieie  punten  heeft  (voor  iedere  wiUdceurige  comUnatk  van  co^  a,  b  en  c  is 
de  fiinctiewaaide  te  beidtoien). 

De  BCONF-ioutine  veilangt  zogdieten  Maitwaarden.  Dit  zijn  waaidor  voor  a,  b,  c  en  coo.  die  als 
b^infNint  voor  het  zodcproces  door  BCONF  ^miikt  worden.  Oindat  BCONF  nia  gamdeeit  dat 
een  i^olwal  opttanum  gevonden  swoidt,  is  het  van  belai^  om  staitwaaiden  le  bqralen  die, 
gesidrstltiieenl  in  ca,  een  zo  kkhi  iiK^lijke  nonn  opteverea 
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62  IMSL-routine  RNUN 


De  RNLIN-routine((^)l  is  speciaal  ontwikkeld  voor  het  oplossen  van  problemen  van  de  volgende 
vorm: 


N-l 

min  I  (yj  -  y(Xi;  a,. 
'=0 


waarbij 

{(Xi,yi)€  R2l0SiSN-l  ) 
gegeven  zijn 

en  a], ....  aj^  onbekende,  te  bepalen  co^fflcienten  voorstellen. 

De  voorwaaide  voor  de  toepassing  van  deze  methode  is  dat  functie  y(.)  bekend  en  condnu  is. 

Het  voordeel  van  deze  methode  is,  dat  het  sooit  problemen  dat  RNLIN  kan  oplossen  oveieenkomt 
met  het  UDB-probleem.  Nadeel  is  dtd  geen  expliciete  intervallen,  waaitwinen  de  coSiTicidnten 
moeten  liggen,  aangegeven  kunnen  worden  Inj  deze  routine.  Ook  RNLIN  heeft  evenals  BCX)NF 
goede  startwaarden  nodig  om  te  kunnen  convergeren  en  ook  RNLIN  garandeett  niet  dat  er  geen 
betere  verzameling  co<!ffici(!nten  bestaat  die  een  kleineie  som  van  de  kwadratische  versdiillen 
gevolg  heeft. 


6J  Keuze  voor  de  optimalizeringsnMthode 

Ofschoon  an  de  voorwaarden  voor  het  gebraik  vn  RNUN  en  BCX>NF  voMan  is,  is  eidge 
vooizichii^ieid  gewenst  De  opiossing  co  als  ftinctie  vn  t,  a,  b,  c  en  <Dd  is  in  het  gdM  random 
b2  >  4  a  c  *  0  wdgedeftniMd.  Er  djn  echter  numerieke  praMemen  te  verwaditen  aai^ezien  de 
ftranules  slecht  geconditioneerd  zijn.  Mtthematisch  gezien  zijn  de  fi>rmules  correct,  maar  ten 
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gevolge  van  het  weiken  met  floating-pmnt  nauwkeurigheid  kunnen  sommige  fiinctie-evaluaUes 
tot  overflow  leiden  (de  noemers  van  mi  en  (03  gaan  naar  nul  als  -  4  a  c  -»  0).  Het  is  deitialve 
aan  te  bevelen  om  in  een  gebiedje  waarvoor  geldt  b2-4ac<e  (e  een  zekere  in  te  steUen 
constante)  de  altijd  welgedefinieerde  (02  te  nemen.  Hieidoor  woidt  editer  kunstmatig  een 
discontinuiteit  ingevoerd.  Robuustheid  is  deitialve  niet  meer  gegiir  ndeeid.  Door  e  zo  klein 
mogelijk  te  kiezen  is  de  kans  dat  de  routine  veikeeide  resultaten  Of^veit  minimaal. 

Een  nadeel  van  RNLIN  is  dat  er  geen  expliciete  grenzen  aangegeven  kunnen  woiden  voor  de 
co€ffici€nten.  RNLIN  maakt  verder  geen  deel  uit  van  de  MATH/IMSL  library,  maar  behoort  tot 
de  STAT/IMSL  library,  die  op  dit  moment  niet  beschikbaar  is.  De  keuze  valt  daaiom  op  BCONF. 
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7  STARTWAARDENBEPALING 


Het  BCONF-algoritme  verlangt  van  de  gebniiker  zogeheten  staitwaaidea  Dit  zijn  waarden  voor 
a,  b.  c  en  (Oq.  die  als  beginpunt  voor  het  zoekproces  van  BCONF  ^bmikt  worden.  Omdat 
BCONF  niet  ganindeert  dat  een  globaal  optimum  gevonden  wordt,  is  bet  zaak  om  startwaaiden  te 
bepalen  die,  gesubstitueeid  in  co,  een  zo  klein  mogelijke  norm  opleveren. 

Voor  de  bepaling  van  een  goede  verzameling  startwaaiden  zijn  er  verschillende  methodcn. 


7.1  Fysische  bepaling  van  de  startwaarden 


Men  kan  een  aantal  vetzamelingen  op  fysische  gionden  gekozen  startwaarden  nemen.  de 
bijbehorende  normen  bepalen  en  de  verzameling  startwaarden,  die  de  kleinste  norm  heeft  als 
invoer  gebruiken  voor  BCONF. 


12  De  kleinste*kwadratenbepaling  van  de  startwaarden 


Deze  methode  gaat  uit  van  de  afgeleiden  van  de  gemeten  hoeksnelheid  naar  de  tijd.  Zij: 


W"'  =  {(tj,  1 0  <  V,  <  ti^  tMax,  i  *  2, 3 . N  - 1;  to  =  0) 


waaibij 

tj  het  tijdsdp  van  de  meting  in  seconden  na  de  eerste  meting, 
co^i'  de  benadeide  eerste  afgeleide  van  de  hodcsnelheid  naar  de  tijd  op  h^  tijdsdp  van 
de  i<^  meting  in  radialen  per  seconde  kwadmat 

Door  de  verzameling  W"'  wordt  een  2^rde  Ueinste-kwadraten  polynoom  gelegd.  Vanwege  de  te 
verwadtten  vorm  van  de  uitrokairve  is  een  benaderir^  van  de  afgeleide  d.m.v.  een  pandnol  een 
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voor  de  hand  liggende  keuze.  De  cofifficiCnien  van  de  parabool  vormen  de  startwaarden  voor  a,  b 
en  c.  De  staitwaarde  voor  (Oq  woidt  gelijk  gekozen  aan  de  eerste  gemeten  hoeksnelheid:  a)™o- 

De  moeten  noodzakelijkerwijs  benadeid  worden,  omdat  de  verandering  in  de  hoeksnelheid 
niet  gemeten  woidt  als  fiinctie  van  de  tijd.  Deze  benadering  voimt  een  nieuw  probleem  en  kan  op 
twee  manieren  worden  aangepakt  die  hieronder  worden  behandeld. 


7  J.1  Lineaire  benadering  van  de  afgeleiden 

Beschouw  de  meeqiuntenverzameling  W™.  De  lineaire  afgeleide  wordt  gedefinieerd  als: 

■  M"i-i 

tom.'*  -  i=1.2 . N-2 

*i+i  •  ^-I 

en 

a)™o’  = 

=w^.2’ 


72.2  Parabolische  benadering  van  de  afgeleiden 

Beschouw  de  meetpuntenverzameling  W™.  De  afgeleiden  worden  benadeid  m.b.v.  de  methodc 
der  kleinste  kwadraten. 

Een  groepje  van  vijf  in  tijd  opeenvolgende  meetpunten  bepaalt  de  afgeleide  in  het  middelste  (3‘) 
meetpunt.  Op  basis  van  de  vijf  meetpunten  woidt  een  kleinste  kwadraten  polynoom  bepaald: 

(12)  «)(t)  =  a  +  Pt  +  yt^  o,p,YeR 

De  afgeleide  van  de  parabool,  met  voor  t  de  tijd  van  het  middelste  meetpunt  gesubstitueeid  is  de 
benaderde  afgeleide  van  het  middelste  punt.  In  foimulevoim: 
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P 

(0™i*  =-  -  ti  i  =  2.  3,...,N-3 

2  Y 


waaibij  P  en  y  de  kleinste-kwadraten  coefficienten  zijn  van  de  parabool  (12),  berekend  m.b.v.  de 
meetpuntendeelveizameling: 

{(tj,  Ciy"j)  I  j  =  i-2,  i-1,  i.  i+1,  i+2) 

ay",'  en  ay"^;^’  worden  lineair  benadeid: 

ay"i^,-oy"i., 

ay"j'=  -  i  =  l.N-2 

•i+l  • 


Voor  de  randwaarden  geldt  tenslotte: 
ay"o'  =  fiy"i' 
ay"N.,'  =<0™N.2‘ 


7J2.3  Methodekeuze  voor  de  benadering  van  de  afgeleiden 

Om  twee  redenen  is  de  parabolische  benadering  van  de  afgeleiden  te  prefereren  boven  de  lineaire 
benadering: 

de  lineaire  benadering  is  gevoeliger  voor  outliers  dan  dc  parabolische  benadering.  Bij  de 
parabolische  benadering  worden  voor  de  benadering  van  de  afgeleide  op  een  zeker  tijdstip 
meer  meetwaaiden  meegenomen  dan  by  de  lineaire  benadering. 
een  parabool  sluit  goed  aan  bij  de  voini  van  de  uitrolcurve. 
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13  Methodekeuze  voor  de  bepaling  van  de  startwaarden 


£>e  bepaling  van  de  staitwaarden  m.b.v.  de  methode  der  kleinste  kwadraten  verdient  de  vooikeur 
boven  de  fysische  bepaling.  De  fysische  bepaling  is  in  zekeie  zin  te  vergelijken  met  bet  prikken 
van  stanwaarden  in  een  4-dimensionale  niimte. 

T.a.v.  de  bepaling  van  de  staitwaaiden  m.b.v.  de  kleinste-kwadraten  methode  kan  opgemeikt 
worden  dat: 

er  piecies  6dn  oplossing  is  voor  a.  ^  en  7  van  (12)  onder  de  aanname  dat  de  curve  niet  te 
benaderen  is  door  een  rechte.  De  aanname  is  noodzakelijk  omdat  anders  A'^A  (zie 
bijlage  B)  singulier  is.  Deze  aanname  is  acceptabel  omdat  uitrolcurves  niet  lineair  zijn. 
de  gevonden  codHlcignten  a,  p  en  y  gegarandeerd  de  kleinste  som  van  de  kwadratische 
verschillen  leveren. 


7.4  Het  fixeren  van  coefTicienten 


Er  kunnen  gevallen  zijn,  waarin  het  wenselijk  is  een  coefficignt  een  vaste  waarde  te  gcvcn  (te 
fixeren).  Indien  bijvoorbeeld  de  c,,-waarde  door  de  fabrikant  van  het  voertuig  aangeleverd  woidt, 
zou  deze  als  een  constante  beschouwd  kunnen  worden.  Misschien  blijkt  het  veider  mogelijk  te 
zijn  op  grond  van  ervaiing  een  vaste  waarde  aan  de  snelheidsonafhankelijke  codfllcient  fj  te 
geven. 

Zowel  de  normfunctie  als  de  bepaling  van  de  startwaarden  moeten  bij  het  fixeren  van 
codfficidnten  aangepast  worden.  De  nonnfunctie  is  eenvoudig  aan  te  passen  door  substihitie  van 
de  vaste  waarde  voor  de  betreffende  parameter  in  de  normfunctie  en  door  het  aanpas.sen  van  het 
aantal  door  de  optimalizeringsmethode  te  bepalen  codfficienten.  De  gevolgen  van  het  fixeren  van 
codflicidnten  voor  de  bepaling  van  de  startwaard^  worden  hieronder  be^rcricen. 
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TS  Gevolgen  fixeren  voor  de  bepaling  van  de  startwaarden 


Indien  geen  van  de  codflici£nten  a,  b  of  c  gefixeerd  woidt,  is  het  toepassen  van  de  medKxle  der 
kleinste  kwadraten  in  feite  niets  anders  dan  het  minimalizeren  van  de  volgende  fiinctie: 

N-l 

(13)  G(a.  b,  c)  =1  (a  tj2  +  b  tj  +  c  -  coP'i)^ 

i=0 

Merk  op  dat  het  fixeren  van  coo  niet  van  invloed  is  op  de  bepaling  van  de  startwaarden  m.b.v.  de 
methode  der  kleinste  kwadraten. 


Indien  a,  b  £n  c  gefixeerd  zijn,  hoeven  er  vanzelfsprekend  ook  geen  staitwaarden  voor  a,  b  en  c 
bepaald  te  worden. 


De  functie  G  is  uit  te  schrijven  als  een  polynoom  van  graad  twee  in  a,  b  en  c: 


(14)  G(a.b,c)=  C^a2 +C,a  + 

C52  b^  +  Cb  b  + 

+CcC  + 

C,b  a  b  +  a  c  +  Cbc  b  c  +  C 


waarbij 


N-l 

C.2  =  I 

i=0 


N*1 

C,  =-2  I 

iaO 


N-l 

Cb2  =  2  tj2 
M) 


N-l 

Cb  =-2  I 

i>0 
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Cc2  =N 

N-1 

Q  =  -  I 

i=0 

N-1 

C.b  =  2  I  ti3 

i=0 

N-1 

C.c  =  2  I 

i=0 

N-1 

Cbc  =  2  X  ti 

i=0 

N-1 

c  =  X 

i=0 


Indien  ddn  of  meer  coefficienten  geflxeenl  zijn,  dan  woiden  hun  waaiden  gesubsUtueerd  in  (13). 
Het  polynoom  dat  oveiblijft  woidt  genoemd. 

Indien  ddn  van  de  codflicidnten  a,  b  of  c  gefixeeid  is,  dan  worden  de  staitwaarden  van  de  niet- 
gerixeerden  bepaald  door  het  stationaiie  punt  van  (een  punt  heet  een  stationair  punt  als  alle 
paittlle  afgeleiden  in  dat  punt  0  zijn).  Dit  stationaire  punt  woidt  bepaald  uit  een  stelsel  van  twee 
lineaiie  vergelijkingen  met  twee  onbekendea  De  coefficientenmatrix  van  het  stelsel  is  mder 
noimale  omstandigheden  niet-singulier.  Immers:  C^»Cu»  Cc2-  Indien  de 
coefTicientenmatrix  singuUer  is,  rest  niets  anders  dan  andere  gefixeeide  waaiden  te  bepalen  en  het 
zoekpioces  te  herhaleru 

Door  de  meetdata  en  gefixeeide  paiameterafhankelijkheid  kan  niet  gegarandeeid  woiden  dat  dit 
stationaire  punt  een  minimum  is.  De  optimalizeringanethode  zal  in  dit  geval  sledite  staitwaaidoi 
aangeboden  krijgen.  Er  zijn  relatief  ingewikkelde  theoriefti  besdiikbaar  om  een  ptdynoom  van 
graad  twee  met  twee  vaiiabelen  te  minimalizeren  biraien  een  zeker  gebied.  Deze  woiden  Mer  itiet 
toegepast,  aangezien  de  coefficiCntenmatrix  zekkn  singulier  zal  zijn. 
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Indien  twee  van  de  drie  coefficignten  a,  b  en  c  gefixeerd  zijn,  blijft  em  polynootn  G,^  van  graad 
twee  over  als  functie  van  de  niet-gefixeerde  parameter.  Door  dit  polynoom  te  optimalizeien 
vinden  we  de  startwaaide  voor  de  niet-gefixeerde  parameter.  De  staitwaaide  voor  de 
niet-gefixeeide  parameter  is  gelijk  aan  die  waaide  waarvoor  het  polynoom  een  optimum 
aanneemt.  Deze  waaide  is  onder  alle  omstandigheden  welgedefinieerd  omdat  C|^  oi 
altijd  ongelijk  aan  nul  zijn.  Er  kan  vanwege  de  afhankelijkheid  van  de  meetdata  en  de  gefixeerde 
co£fficii!nten  niet  gegarandeeid  woiden  dat  dit  optimum  altijd  een  minimum  is.  Als  bet  optimum 
geen  minimum  is,  zal  de  optimalizeringsmethode  slechte  staitwaaiden  aangeboden  krijgen. 
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8  CONVERGENTIE  EN  UNICITEIT 


Zij  f  de  normfiinctie: 

f  :  R4  R 

(a,  b,  c,  (Oq)  -*  f(a,  b,  c,  Mq) 

en  zij  ( a'^.  b^,  c'^,  <Oo'^)  ccn  lokaal  minimum  van  f. 

T.a.v.  de  nauwkeurigheid  en  de  uniciteit  kunnen  de  volgende  vragen  gesteld  woiden: 

1.  Heeft  f  een  oplossing  met  de  codfficignten  binnen  de  opgegeven  intervallen  ?  Zo  ja,  is  er 
dan  precies  66n  oplossing  ? 

2.  Leidt  de  optimalizeringsmethode  tot  een  oplossing,  zo  mogelijk  tot  d6  oplossing  ? 

3.  Als  (a^,  b^,  c\  (Oo^)  bepaald  is,  tot  op  hoevecl  decimalen  nauwkeurig  zijn  de 
componenten  bepaald  ? 

ad  1.  Omdat  f  een  in  bet  zoekgebied  welgedefinieeide,  kwadratische  functie  is,  is  er  een 
oplossing  (f  2  0 1 V  a,  b,  c,  (Oo  e  [zoekgebied])  binnen  het  zoekgebied.  Niet  gegarandeerd 
is  dat  er  slechts  6dn  oplossing  is. 

In  het  geval  f  een  polynoom  is  dan  zijn  er  methoden  om  aan  te  tonen  dat  f  binnen  een 
interval  ddn  of  meer  extremen  heeft.  E6n  van  de  methoden  (SturmK^)])  weikt  mtt  behu^) 
van  differentiatie  'zolang  als  nodig'.  Het  polynoom  wordt  'opgedifTerentieenl'  en  ieder  van 
de  afgelekten  wordt  afzondeilijk  ondenocht  op  extremen  en  nulpunten.  In  66n  dimensie 
kan  dit  al  veel  work  inhouden. 

Bij  toepassing  van  DBCONF  met  de  of^ossing  to  verweikt  in  f  is  de  bovengaioemde 
methode  niet  toqMsbaar.  f  is  immers  geen  polynoom.  Mij  zijn,  behalve  de  gewone 
analytisdie  methoden  (ftmctie-onderzoek),  geen  methoden  bekend  om  een  deigeU^ 
4-dimensionaal  pioUeem  am  te  piddtan.  Ocric  zonder  gebraflt  vm  oplossta^  co  te  het 
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bepalen  van  de  aanwezi^id  van  66n  of  meer  extiemen  een  lasUge  zaak.  De  fiinctie  f  is 
inuners  afhankelijk  van  meetwaaiden,  terwijl  een  van  de  meetwaarden-(mafhaidttli^ 
uitspraak  wenselijk  is.  Schnabeli<^>i  stelt  over  dit  type  probkem  dat  'men  veelal  genoegen 
moet  nemen  met  een  door  de  computer  aangebodoi  oplossingf.  De  ojdossing  dient  dan 
wel  vergezeld  te  gaan  met  een  norm  of  maat,  zodat  op  grand  daarvan  gezegd  kan  woiden 
of  men  wel  of  niet  tevieden  is  met  bet  lesultaat 

ad  2.  Omdat  de  startwaaiden  binnen  de  opgegeven  intervallen  moeten  liggoi  vindt  de 
optimalizeringsmethode  altijd  een  ofdossing.  Het  gebniik  van  BCONF  garandeett  niet  dat 
bet  beste  minimum  gevonden  woidL 

ad  3.  Om  een  uitspraak  over  het  aantal  goede  decimalen  te  doen  is  ondeizodc  van  de  paitifle 
afgeleiden  noodzakelijk.  De  eeiste  oide  paitieie  afgeleiden  vertellen  On  de  eerste  oide) 
wat  de  verandering  in  de  iunctiewaaide  is  indien  een  parameter  een  klein  bcetje  varieeit 
Om  tot  een  algemene  uitspraak  te  komen  die  zich  niet  beperkt  tot  een  spedfiek  punt 
(a,  b,  c,  coo),  is  het  noodzakelijk  om  de  eerste  oide  codfRcienten  af  te  schatten.  De  partible 
afgeleiden  van  de  oplossingen  zijn  in  dit  specifieke  geval  helaas  onweikbaar  gioot. 

Er  kan  geconcludeerd  woiden  dat  er  altijd  een  oplossing  gevonden  woidt;  er  valt  echier  weinig  te 

zeggen  over  de  nauwkeurigheid  en  de  uniciteit  van  de  oplossing  (zie  hoofdstuk  10  voor  de  testen). 
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9  HEX  SCHATTERALGORITME 


De  werking  van  het  schatteralgoritme  in  het  UDB  project  kan  als  volgt  samengevat  woiden: 

1 .  Bepaling  van  de  afgeleiden  met  de  methode  der  kleinste  kwadraten. 

2.  Bepaling  van  de  staitwaaiden  m.b.v.  de  methode  der  kleinste  kwadraten.  Startwaaiden 
mogen  gefixeeid  worden. 

3.  Bepaling  van  de  optimale  colifricienten  m.b.v.  de  IMSL-ioutine  DBCONF,  gebruik 
makoide  van  de  exacte  oplossingen,  waaibij  (Dq  een  parameter  is. 


Voor  de  implementatie  van  de  schatteralgoritme  kunnen  de  volgende  punten  opgemeikt  worden: 

De  startwaaiden  hoeven  niet  in  door  de  gebruiker  aangegeven  intervallen  te  liggea  Indien 
een  startwaarde  niet  in  het  interval  voor  die  specifieke  coefficient  ligt,  dan  woidt  deze 
gelijk  gestdd  aan  de  dichtstbijzijnde  intervalgrens. 

Het  oplossen  van  de  stelsels  lineaire  vergelijkingen  Lb.v.  de  bepaling  van  de  afgddden 
van  de  hoeksnelheid  behoeft  geen  pioblemen  op  te  leveren  aangezien  de  tijd  tussen  twee 
opeenvolgende  metingen  voldoende  lang  is  om  singulariteit  van  A’^A  uit  te  sluiiea 
Omdat  BCX)NF  yfni're  d^erences  gebruikt  om  de  partieie  afgeleiden  van  de  normfunctie 
te  benadeitn  is  het  wenselijk  de  dubbele  precisie  versie  van  BCX)NF  te  gdxuiken: 
DBCXDNF. 

T.b.v.  de  precisie  van  de  lesulUtten  is  het  wenselijk  DBCONF  factorcn  te  laien  beiekenea 
waaimee  de  startwaaiden  vermenigvuldigd  moeten  worden  om  de  optimak 
mathematische  coCRicienten  te  berekenen.  Dit  vooikomt  dat  DBCONF  mtt  heel  grote  rtf 
juist  hed  kldne  coCfficienten  lekent,  waardoor  afhXNlfouten  een  groteie  rol  kunnen 
spdea 

Er  dimien  voorzorgen  genomen  te  worden  om  overflow  van  m  te  vooikomca 
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10  TESTEN  SCHATTERALGORITME 


10.1  Het  dod  van  de  testen 

Een  van  de  problemen  bij  de  bepaling  van  de  optimale  waanlen  van  de  parameters  fi.  oi  is. 
dat  de  te  minimalizeren  fiinctie  niet  noodzakelijkerwijs  slechts  ddn  minimum  in  het  zodcgebied 
behoeft  te  hebben.  Het  is  derhalve  niet  uitgesloten  dat  het  zoekalgoiitme  een  veikcerd  (met 
optimaal)  minimum  vindt,  hetgeen  onjuiste  coefficiiinten  tot  gevolg  heeft.  Het  doel  van  de  testen 
is  aan  te  tonen  dat  het  schatteralgoritme  m.b.v.  de  door  het  schatteralgoritme  berekende 
startwaatden  in  praktisch  aUe  gevallen  het  beste  optimum  vindt 


10.2  Testprocedure 

Aangezien  het  praktisch  onmogelijk  is  op  analytische  wijze  te  bewijzen  dat  er  geen  andere  bctere 
(kleinere)  minima  bestaan,  is  gepoogd  op  munerieke  wijze  een  indnik  te  krijgen  van  het  karakter 
van  de  noimfiinctie. 

Daaitoe  woiden  met  het  schatteralgoritme  de  optimale  coCflicienten  bepaald  behorende  bij  een 
uitrolcurve  uit  de  praktijk.  Deze  c^mak  coSfficiCnten  leggen  een  collectie  startwaarden  vast: 

{ (a  *  f,.  p  •  fj.  Y *  c,)  I  a.  p.  Y»  0.1. 1.0. 10.0 ) 

Op  ieder  tupel  startwaarden  (f],  fj,  c^)  is  nu  het  sdutteralgoritme  toegqwst.  waaibij  de 
stattwaarde  voor  aio  steeds  gelijk  gesteld  wordt  aan  de  eerste  gemeten  hoeksneDiekl.  De 
berekeidng  van  de  startwaarden  door  het  achatteralgotitme  wordt  in  dit  geval  niet  uitgevoerd. 
Door  de  omvmg  van  de  ooHecde  (27  tqtels  startwaarden)  kan  aat^enomen  wmden  (tat  de 
GoOfflcientenvenamelinf  met  de  kleinste  norm  het  Abaoluie  mininHim  voimt  bbmen  bet 
zoekgrtiied.  Als  deze  coeffkWnienventmneMng  ovoeenkomt  mm  de  door  het  acbtatendgorilme 


i 

I 
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(m^  de  startwaaidenbepaling)  berekende  opdmale  coefficienten.  dan  heeft  het  sdiatteialgoritme 
het  gewenste  absolute  minimum  blijkbaar  bepaald  voor  de  specifieke  uitiolcurve. 

Deze  testpiocedure  woidt  uitgevoerd  met  verschillende  uitiolcufves.  Het  volstaat  ni^  om  aan  te 
tonen  dat  in  een  speciiiek  geval  het  schatteralgoritme  het  absolute  minimum  vindt 


lOJ  Testdata 


Het  schatteralgoritme  is  getest  met  een  vijftal  uitrolcurves  uit  de  praktijk: 

UitrollA 

UitiolIB 

Uitiol2A 

UitiDl2B 

Land 


De  laatste  uitiolcurve  is  afkomstig  van  een  Landiover.  Alle  andere  zijn  gemaakt  met  de  DAF  YA- 
4440.  De  uitrolcurves  lA  en  2A  zijn  in  dezelfde  (wind-)  richting  uitgevoerd.  De  richting  waarin 
IB  en  2B  uitgevoerd  zijn  was  tegengesteld  aan  die  van  1 A  en  2A.  De  windsnelheid  is  niet  bckend. 

Hieronder  volgen  de  door  het  schatteralgoritme  gebruikte  technische  gegevens  van  de  voeituigen: 


Voertuig 

Massa[kg] 

D]mamisdie  wielradius  [m] 
Effectieve  oppervldc  [m^] 


YA-4440  Landiover 

6480  1750 

0J25  0.330 

4.1625  2.3200 


Het  zoekgdiied  votwoptimale  coeflicUlnten  was  voor  alle  uhrolcurves: 
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0.0001  ^  f,  ^  10.0000 
0.000001  10.0000 

0.00001  10.0000 

10.0  Vo  ^100.0  km/h 


10.4  Toelichting  beschryving  testresultaten 


Van  elke  uitrolcurve  is  een  grafiek  met  de  snelheid  als  hmctie  van  de  tijd  6n  een  grafiek  met  het 
snelheidsverschil  tussen  de  berekende  uitrolcurve  en  de  gemeten  uitrolcurve  gemaakt  Een  tabel 
met  tupels  en  bijbehorende  optimale  co€f!ici6nten  is  gegeven  en  de  test  is  van  commentaar 
vooTzien. 

De  graOeken  geven  een  visuele  indruk  van  de  kwaliteit  van  de  fit.  Indien  de  berekende  en  de 
gemeten  uitrolcurve  overeenkomen  dan  is  er  sprake  van  een  goede  fit  De  coOrdinaten  van  de 
berekende  uitrolcurve  zijn  verktegen  door  de  optimale  coefficiSnten  te  substitueren  in  de  exacte 
opiossing,  waardoor  een  fiiiKtie  van  aileen  de  tijd  ontstaat. 
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Verklaring  legends  uitrolcurve-grafiekoi  (zie  b.v.  iiguur  1  op  pagina  40) 


file.txt 

de  naam  van  de  file  van  de  uitrolcurve 

massa 

de  massa  van  bet  voertuig  [kg] 

DWR 

de  dynamische  wielradius  [m] 

EffA 

de  effectieve  oppervlakte  [m^j 

fi 

de  optimale  coefficient  f] 

f2 

de  optimale  coefficient  f2 

de  optimale  coefficient  c^ 

Wo 

de  optimale  coefficient  (Do 

Delta 

de  bijbehorende  A  (»  noim) 

Origineel 

De  lijn  'Origineel*  is  de  polylijn  van  de  gemeten  snelheid  als  fimctie  van  de  tijd 

Resultaat 

De  lijn  'Resultaat*  is  de  polylijn  van  de  berekende  snelheid  als  functie  van  de  tijd. 

&  is  gedeOnieerd  door  (zie  (1 1)); 
A  =  V(F(a.b.c,an,)/N) 
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Verklaring  van  de  kolommen  van  de  tabel  (zie  b.v.  label  Uitroll  A  op  pagina  41): 


n  Identilicatie  nummer  van  de  tupel 
fi  Staitwaarde  van  f] 
f2  Staitwaaide  van  £2 
Staitwaarde  van  c,, 

oiq  Hoeksnellieid  op  het  tijdstip  van  de  eerste  meting  [rad  /  s] 

A  De  norm  behorende  bij  de  staitwaarden  [rad  /  s].  Een  zeer  hoge  A  geeft  aan  dat  de 
staitwaarden  zeer  slecht  zijn.  Het  is  derhalve  onwaarschijnlijk  dat  dergelijke 
staitwaarden  berekend  zullen  woiden  door  de  methode  der  kleinste  kwadratea 
Opt_fi  Optimale  waaide  ft 
Opt_f2  Optimale  waarde  f2 
Opt_Cv,  Optimale  waarde 

Opt_(i]()  Optimale  hoeksnelheid  op  het  tijdstip  van  de  eerste  meting  [rad  /  s] 

Opt_A  Optimale  A  behorende  bij  de  optimale  waaiden  [rad  /  s] 

WC  Warning  Code.  Een  van  DBCONF  afkomstige  waarschuwingscode. 

0  geen  problemen  tijdens  zoekproces 
3  maximum  aantal  iteraties  overschieden. 

8  mogelijk  niet  optimale  minimum  gevonden 

Andeie  waarden  zijn  00k  mogelijk,  maar  hebben  zich  niet  voorgedaan  tijdens  de 
testen. 

Iter  Het  aantal  iteraties.  DBCONF  is  een  iteratief  proces.  Het  ingestelde  maximum  aantal 
iteraties  bedraagt  100. 


De  kolommen  voor  tag  zijn  toegevoegd  aan  de  tabel  aangezioi  dit  tijdens  het 
opUmalizeringspioces  een  gewone  parameter  is.  Een  extreme  waarde  voor  (Oo  kan  eat  verklaring 
vormen  voor  een  grote  A. 

In  de  tabel  is  Opt_A  de  meest  interessante  kolom.  Deze  kolom  kan  het  aannemelijk  makra  dat  de 
te  minimalizeien  fimctie  precies  ddn  minimum  heeft  oi  dus  dat  de  param^ers  f],  f2  en  c«  uniek 
bepaald  zijn.  Daar  waar  Opt_A  alwijkt  v«i  hetgeoi  tq)  basis  van  de  beste  waaiden  voor  de 
parameters  f],  f2  en  verwadn  had  mogen  woiden,  woidt  een  verklaring  gegeven.  Om 
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onderscheid  te  maken  tussen  makkelijke  en  moeilijke  gevallen  zijn  de  tupels  geclassificeerd.  Er 
zijn  3  klassen:  goede.  madge  en  sledite  tupels. 

Een  tupel  is  goed  indien  de  by  de  opdmale  parameteis  behoiende  A  de  kleinste  is  van  de  hele 
kolom. 

Een  tupel  is  matig  indien  de  bij  de  optimale  parameters  behoiende  A  'niet  veel'  afwijkt  (maximaal 
een  factor  2)  van  de  kleinste  A  uit  de  kolom. 

Een  tupel  is  slecht  indien  het  maximum  aantal  iteraties  uitgevoerd  was  of  als  A  'groot'  is  (d.i.: 
meer  dan  een  factor  2  afwijkt  van  de  kleinste  A). 

Deze  classificatie  is  ook  tenig  te  vinden  in  de  tabel: 

De  slechte  tupels  zijn  voorzien  van  een  ster 
De  Opt_A  van  de  matige  tupels  is  onderstreept. 

De  goede  tupels  hebben  geen  ster  en  ook  is  de  Opt_A  niet  onderstreept. 


De  matige  tupels  zouden  bij  ontbreken  van  andere  resultaten  kunnen  leiden  tot  incorrecte 
resultaten. 


Indieu  blijkt  dat  er  voor  de  kleinste  A's  slechts  €6n  unieke  optimale  parameterveizameliitg  bestaat, 
dan  is  het  aannemelijk  dat  de  te  minimalizeren  functie  slechts  6^n  minimum  heeft.  Blijkt  echter 
dat  er  meer  niet-identieke  parameterveizamelingen  bestaan  voor  de  kleinste  A's,  dan  is  het 
aannemelijk  dat  de  te  minimalizeren  fimetie  verschillende  minima  heeft.  De  waaiden  voor  f^,  (2 
en  Cw  zijn  dan  afhankelijk  van  de  beginschatting  en  niet  uniek  bepaald. 


Bij  het  commentaar  wonft  ^cifiek  aandacht  gegeven  aan: 

De  ligging  van  de  origineel-  en  resultaatpolylijn  t.o.v.  elkaar. 

De  veihouding  Opt_A  /  Opt_(ao.  Des  te  kleiner,  des  te  meer  lijkt  de  resultaatcurve  op  de 
origineelcurve. 

Beschouwing  vvi  de  slechte  en  matige  tupds. 

Unicitdt  miidmun. 

Het  veigdijk  van  het  minimum  bepaald  door  het  schatteralgoritme  (met 
staitwaanlenbepaling)  en  het  bene  minimum  uit  de  tabd. 
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Tabel  1.  UitrollA  coeffidenten  bepaald  uitgaande  van  verschillende  startwaarden 


n 

«1 

^2 

•o 

A 

Opt^fl 

Opt^fj 

Opt_c^ 

Opt_«o 

Opt_A 

we 

lt«r 

1 

0.00076 

0.00003 

0.04300 

46.270 

26.9211  1 

0.00076 

0.00001 

0.00001 

38.76? 

2.5736 

3 

100* 

2 

0.00076 

0.00003 

0.43000 

46.270 

11,4670  1 

0 . 00760 

0.00031 

D.435S5 

46.358 

0.2066 

0 

23 

3 

0.00076 

0.00003 

4.30000 

46.270 

10.8076  1 

0.00760 

0.00031 

0.43555 

46.358 

0.2066 

0 

31 

4 

0.00076 

0.00030 

0.04300 

46.270 

17.4846  1 

0.00010 

0.00002 

0.00008 

38.884 

2.5448 

3 

200* 

5 

0.00076 

0.00030 

0.4JOOO 

46.270 

7.4696  1 

0.00760 

0.00031 

0.43555 

46.358 

0.2066 

0 

19 

6 

0.00076 

0.00030 

4.30000 

46.270 

11.3418  1 

0.00760 

0.00031 

0.43555 

46.358 

0.2066 

0 

35 

7 

0.00076 

0.00300 

0.04300 

46.270 

9.41S2  1 

0.00760 

0.00031 

0.43555 

46.358 

0.2066 

0 

28 

8 

0.00076 

0.00300 

0.43000 

46.270 

10.S161  1 

0.00760 

0.00031 

0.43555 

46.358 

0.2066 

0 

32 

9 

0.00076 

0.00300 

4.30000 

46.270 

15.0816  t 

0.00760 

0.00031 

0.43555 

46.358 

0.2066 

0 

47 

10 

0.00760 

0.00003 

0.04300 

46.270 

14.8871  1 

0.00760 

0.00031 

0.43555 

46.358 

0.2066 

0 

28 

11 

0.00760 

0.00003 

0.43000 

46.270 

3.7762  1 

0.00760 

0.00031 

0.43555 

46.358 

0.2066 

0 

15 

12 

0.00760 

0.00003 

4.30000 

46.270 

16.3080  1 

0.00760 

0.00031 

0.43555 

46.358 

0.2066 

0 

27 

13 

0.00760 

0.00030 

0.04300 

46.270 

8.0358  1 

0.00760 

0.00031 

0.43555 

46.358 

0.2066 

0 

26 

14 

0.00760 

0.00030 

0.43000 

46.270 

0.2398  1 

0.00760 

0.00031 

0.43555 

46.351 

0.2066 

0 

9 

13 

0-00760 

0.00030 

4.30000 

46.270 

16.7905  1 

0.00760 

0.00031 

0.43555 

46.351 

0.2066 

0 

33 

14 

0.00760 

0.00300 

0.04300 

46.270 

12.6080  1 

0.00760 

0.00031 

0.43555 

46.351 

0.2066 

0 

27 

17 

0.00760 

0.00300 

0.43000 

46.270 

13.5161  1 

0.00760 

0.00031 

0.43555 

46.358 

0.2066 

0 

24 

11 

0.00760 

0.00300 

4.30000 

46.270 

17.7035  1 

0.00760 

0.00031 

0.43555 

46.358 

0.2066 

0 

34 

19 

0.07600 

0.00003 

0.04300 

46.270 

114.3363  1 

0.00010 

0.00009 

0.00176 

39.321 

2.3741 

3 

100* 

20 

0.07600 

0.00003 

0.43000 

46.270 

1070.3910  1 

0.00010 

0.00001 

0.00001 

38.763 

2.5767 

3 

100* 

21 

0.07600 

0.00003 

4.30000 

46.270 

1102.7655  1 

0.2091$ 

0.00041 

10.00000 

10.000 

176. 1522 

1 

25* 

22 

0.07600 

0.00030 

0.04300 

46.270 

101.6421  1 

0.00010 

0.00001 

0.00003 

38.770 

2.5665 

3 

100* 

23 

0.07600 

0.00030 

0.43000 

46.270 

1220.2246  ( 

0.00010 

0.00002 

0.00002 

38.772 

2.5706 

3 

100* 

24 

0.07600 

0.00030 

4.30000 

46.270 

366.8917  1 

0.00760 

0.00031 

0.43555 

46.358 

0.2066 

0 

48 

25 

0.07600 

0.00300 

0.04300 

46.270 

48.5530  1 

0.00760 

0.00031 

0.43555 

46.358 

0.2066 

0 

33 

26 

0.07600 

0.00300 

0.43000 

46.270 

57.8503  I 

0.00760 

0.00031 

0.43555 

46.358 

0.2066 

0 

33 

27 

0.07600 

0.00300 

4.30000 

46.270 

269.1175  1 

0.09546 

0.00128 

10.00000 

28.846 

145.2587 

3 

100* 

Oi^neilungen; 

De  grafieken  tonen  dat  de  origineel-  en  resultaatpolylijn  keurig  op  elkaar  liggen. 

Voor  de  kleinste  Opt_A  bedraagt  de  veihouding  Opt_A  /  Opt_<Do  (xigeveer.  0.5  %. 
(0.2066  /  46.385) 

Alle  slechte  tupels  onderscheiden  zich  duidelijk  door  de  sterk  afwijkende  Opt_A.  Gevolg: 
er  resteren  enkel  goede  tupels  en  er  is  uniek  minimum  gevondea 
Het  schatteralgoritme  (met  startwaaidenbepaling)  heeft  het  absolute  minimum  in  bet 
zoekgebied  gevondea 
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Tabel  2.  UitrollB  coefficienten  bepaald  uitgaande  van  verschillende  startwaarden 


n 

'l 

^2 

Si 

•o 

A 

Opt_fi 

Opt^fj 

Opt_C^ 

Opt_«0 

Opt_A 

NC 

zt«r 

X 

0.00076 

0.00001 

0.07070 

46.270 

24.6936 

0.00010 

0.00001 

0.00001 

37.643 

2.6625 

3 

100* 

2 

0.00076 

0.00001 

0.70700 

46.270 

7.0529 

0.007S9 

0.00000 

0.70642 

46.762 

0.5944 

0 

15 

3 

0.00076 

0.00001 

7.07000 

46.270 

13.0664 

0.00742 

0.00000 

0.72797 

47.317 

0.6139 

0 

37 

4 

0.00076 

0.00010 

0.07070 

46.270 

22.4316 

0.00010 

0.00002 

0.00022 

37.610 

2.6460 

3 

100* 

i 

0.00076 

0.00010 

0.70700 

46.270 

6.0393 

0.00759 

0.00000 

0.70642 

46.762 

0.5944 

0 

13 

C 

0.00076 

0.00010 

7.07000 

46.270 

13.1667 

0.00759 

0.00000 

0.70642 

46.762 

0.5944 

0 

40 

7 

0.00076 

0.00100 

0.07070 

46.270 

3.6006 

0.00759 

0.00000 

0,70642 

46.762 

0.5944 

0 

19 

t 

0.00076 

0.00100 

0.70700 

46.270 

3.4232 

0.00759 

0.00000 

0.70642 

46.762 

0.5944 

0 

16 

9 

0.00076 

0.00100 

7.07000 

46.270 

14.2477 

0.00759 

0.00000 

0.70642 

46.762 

0.5944 

0 

36 

10 

0.00760 

0.00001 

0.07070 

46.270 

14.0662 

0.00793 

0.00000 

0.63211 

45.667 

0.6951 

0 

16 

u 

0.00760 

0.00001 

0.70700 

46.270 

0.6637 

0.00759 

0,00000 

0.70642 

46.762 

0.5944 

0 

9 

12 

0.00760 

0.00001 

7.07000 

46.270 

67.2049 

0.01135 

0.00000 

10.00000 

65.746 

42.0363 

0 

24* 

13 

0.00760 

0.00010 

0.07070 

46.270 

11.6430 

0.00759 

0.00000 

0.70642 

46.762 

0.5944 

0 

24 

14 

0.00760 

0.00010 

0.70700 

46.270 

1.3694 

0.00759 

0.00000 

0.70642 

46.762 

0.5944 

0 

6 

13 

0.00760 

0.00010 

7.07000 

46.270 

71.6363 

0.01135 

0.00000 

10.00000 

65.746 

42.0363 

0 

24* 

16 

0.00760 

0.00100 

0.07070 

46.270 

3.5763 

0.00759 

0.00000 

0.70642 

46.762 

0.5944 

0 

26 

17 

0.00760 

0.00100 

0.70700 

46.270 

6.1330 

0.00759 

0.00000 

0.70642 

46.762 

0.5944 

0 

12 

11 

0.00760 

O.QQIOO 

7.07000 

46.270 

195.0702 

0.01069 

0.00000 

10.00000 

65.765 

42.0116 

0 

36* 

19 

0.07600 

0.00001 

0.07070 

46.270 

131.6747 

0.00010 

0.00001 

0.00002 

37.661 

2.6766 

3 

100* 

20 

0.07600 

0.00001 

0.70700 

46.270 

1167,3297 

0.47166 

0.00000 

0.00007 

10.000 

14.2296 

6 

15* 

21 

0.07600 

0.00001 

7.07000 

46.270 

1131.3095 

0.00755 

0.00000 

0.70620 

46.765 

0.5966 

0 

41 

22 

0.07600 

0.00010 

0.07070 

46.270 

125.4606 

0.00753 

0.00000 

0.66956 

46.245 

0.6254 

0 

29 

23 

0.07600 

O.OOOlO 

0.70700 

46.270 

547.3117 

0.00777 

0.00000 

0.66134 

46.467 

0.6030 

0 

32 

24 

0.07600 

0.00010 

7.07000 

46,270 

794.3610 

0.00049 

0.00016 

1.32750 

49.223 

3.1717 

3 

100* 

23 

0. 07600 

0.00100 

0.07070 

46.270 

63,1699 

0.00759 

0.00000 

0.70642 

46.762 

0.5944 

0 

31 

26 

0.07600 

0.00100 

0.70700 

46.270 

666.1066 

0.00759 

0.00000 

0.70642 

46.762 

0.5944 

0 

43 

27 

0.07600 

0.00100 

7.07000 

46.270 

137.4016 

0.00759 

0.00000 

0.70642 

46.762 

0.5944 

0 

42 

Conunoittaar. 

De  graTieken  tonen  dat  de  origineel-  en  6c  resultaaqwlylijn  relatief  veel  versdiilloi  rond 
tijdstip  118  seconden.  Navraag  leeide  dat  de  afwijking  wellicht  een  gevolg  is  van 
(flauwe)  boditen  aan  het  eind  van  de  testweg.  Deze  oefenen  een  extra  kiacht  uit  op  het 
voeituig.  De  omclusie  die  hieruit  getrdcken  mag  woiden  is,  dat  de  berdcende  optimak 
codfRciftiten  niet  bdrouwtNur  zijn  doo*  de  slechte  originedpdylijn. 

Voor  de  kleinste  Opt_A  bedraagt  de  veihouding  C)pt_A/Opt_(ioo  ongeveen  1.3  %. 
(0.5944/46.762) 

AUe  dechte  tupels  ondeischeiden  zich  duiddijk  door  de  steik  afWijkende  Opt_A. 
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De  matige  tupds  kunnen  niet  afgewezen  worden.  Het  versdiil  tussen  de  Opt_A's  van  de 
matige  tupds  is  erg  klein:  ongeveer  S%  (0.03/0.60).  De  uitspraak  dat  hietmee 
aangetoond  is  dat  de  optimale  coSffidenta)  staitwaaide-afiiankelijk  zoudoi  zijn  is  editer 
niet  te  maken  aangezien  de  origineelpolylijn  slecht  is.  Bovendien  liggen  de  bij  de  matige 
en  goede  A's  behoiende  optimale  coeffidditen  niet  ver  van  dkaar  af.  Groteie 
tekenpiecisie  of  een  kleiner  stopcriterium  voor  DBCONF  zouden  mogelijk  tot  gevolg 
hebben  dat  convergentie  plaats  vindt  naar  ddn  minimum,  als  er  in  weikelijkheid  lond  het 
absolute  minimum  een  gebiedje  is  waar  de  punten  byna  stadonair  zijn. 

Een  fysische  reden  om  twijfelen  aan  de  origineelpolylijn  is  de  oniealistische  waaide  voor 
Opt.fj;  0.00000. 

Ondanks  de  slechte  origineelpolylijn  komt  het  minimum  geleverd  door  het 
schatteralgoritme  (met  startwaardenbepaling)  oveieen  met  het  absolute  minimum  in  het 
zoekgebied  volgens  de  tabel. 


i 
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Tabel  3.  UitroI2A  co^ffidenten  b^aaid  uitgaande  van  vo'scbillende  startwaarden 


n 

fl 

^2 

*0 

A 

Opt^e^ 

opt_«o 

OptjA 

we  lt*r 

1 

0.00077 

0.00003 

0.04200 

46.063 

26.9102 

i 

0.00077 

0.00001 

0.00001 

36 

590 

2.5036 

3 

100* 

2 

0.00077 

0.00003 

0.42000 

46.063 

11.7946 

1 

0.00773 

0.00031 

0.42722 

45 

930 

0.1630 

0 

23 

3 

0.00077 

0.00003 

4.20000 

46.063 

10.5323 

1 

0.00773 

0.00031 

0.42722 

45 

930 

0.1630 

0 

26 

4 

0.00077 

0.00030 

0.04200 

46.063 

17.5525 

t 

0.06010 

0.00002 

0.00010 

36 

702 

2.4724 

3 

100* 

5 

0.00077 

0.00030 

0.42000 

46.063 

7.7401 

1 

0.00773 

0.00031 

0.42722 

45 

930 

0.1630 

D 

20 

€ 

0.00077 

0.00030 

4.20000 

46.063 

11.0753 

1 

0.00773 

0.00031 

0.42722 

45 

930 

0.1630 

0 

36 

1 

0.00077 

0.00300 

0.04200 

46.063 

9.2769 

I 

0.00773 

0.00031 

0.42722 

45 

930 

0.1630 

0 

29 

• 

0.00077 

0.00300 

0.42000 

46.063 

10.3564 

1 

0.00773 

0.00031 

0.42722 

45 

930 

0.1630 

0 

35 

9 

0.00077 

0.00300 

4.20000 

46.063 

14.6961 

1 

0.00773 

0.00031 

0.42722 

45 

930 

0.1630 

0 

42 

10 

0.00770 

0.00003 

0.04200 

46.063 

14.6261 

1 

0.00773 

0.00031 

0.42722 

45 

930 

0.1630 

0 

27 

u 

0.00770 

0.00003 

0.42000 

46.063 

4.0079 

1 

0.00773 

0.00031 

0.42722 

45 

930 

0.1630 

0 

17 

12 

0.00770 

0.00003 

4.20000 

46.063 

15.7666 

1 

0.00773 

0.00031 

0.42722 

45 

930 

0.1630 

0 

26 

13 

0.00770 

0.00030 

0.04200 

46.063 

6.0533 

1 

0.00773 

0.00031 

0.42722 

45 

930 

0.2630 

0 

24 

14 

0.00770 

0.00030 

0.42000 

46.063 

0.3764 

1 

0.00773 

0.00031 

0.42722 

45 

930 

0.1630 

0 

10 

15 

0.00770 

0.00030 

4.20000 

46.063 

16.2696 

1 

0.00773 

0.00031 

0.42722 

45 

930 

0.1630 

0 

33 

16 

0.00770 

0.00300 

0.04200 

46.063 

12.5055 

I 

0.00773 

0.00031 

0.42722 

45 

930 

0.1630 

0 

26 

17 

0.00770 

0.00300 

0.42000 

46.063 

13.3952 

1 

0.00773 

0.00031 

0.42722 

45 

930 

0.1630 

0 

24 

It 

0.00770 

0.00300 

4.20000 

46.063 

17.5407 

1 

0.00773 

0.00031 

0.42722 

45 

930 

0.1630 

0 

29 

If 

0.07700 

0.00003 

0.04200 

46.063 

114.2764 

1 

0.00010 

0.00009 

0.00200 

39 

260 

2.3074 

3 

100* 

20 

0.07700 

0.00003 

0.42000 

46.063 

1212.7660 

1 

0.00010 

0.00001 

0.00001 

36 

565 

2.5071 

3 

100* 

21 

0.07700 

0.00003 

4.20000 

46.063 

246.2373 

1 

0.01264 

0.00000 

10.00000 

10 

000 

46.4252 

6 

41* 

22 

0.07700 

0.00030 

0.04200 

46.063 

101.7461 

1 

0.00010 

0.00001 

0.00002 

36 

606 

2.4993 

3 

100* 

23 

0,07700 

0.00030 

0.42000 

46.063 

1233.6154 

r 

0.00010 

0.00026 

1 . 92796 

13 

526 

14.5420 

0 

5* 

24 

0.07700 

0.00030 

4.20000 

46.063 

364.2105 

1 

0.22146 

0.00000 

10.00000 

13 

365 

107.5965 

0 

16* 

25 

0.07700 

0.00300 

0.04200 

46.063 

46.6409 

1 

0.00773 

0.00031 

0.42722 

45 

930 

0.1630 

0 

35 

26 

0.07700 

0.00300 

0.42000 

46.063 

57.9394 

1 

0.00691 

0.00000 

0.61695 

46 

319 

0.2279 

6 

33 

27 

0.07700 

0.00300 

4.20000 

46.063 

625.9624 

1 

0.20909 

c.coooo 

10.00000 

34 

063 

172.6129 

0 

43 

Opmeikingen: 

De  graTieken  umen  dat  de  origineel-  en  resultaatpolylijn  keurig  op  elkaar  liggen. 

Voor  de  kldnste  Opt_A  bedraagt  de  veifiouding  Opi_A/Optjaio  ongeveen  0.4  %. 
(0.1630  /46.358) 

Alle  slechte  tupels  onderscheiden  zkh  duidelijk  door  de  steik  afwijkende  Opt_A. 

De  enige  matige  tupel  (26)  kan  op  fysische  grondoi  uitgesloten  wmden  aangezien  OpL.f2 
een  findwaarde  van  bet  zodcgebied  aanneemt 

Gevolg:  er  resieren  enkd  goede  tupeb  en  O’ is  66n  unidc  minimum  gevondoL 

Het  schaaeialgoritme  (met  staitwaaixknbepaling)  heefk  bet  abaolute  minimum  in  bet 

zoekgebied  gevondea 


UitroUB  test 


UITR0L2B.TXT 

Massa  MO 


i 
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Tabel  4.  Uitrol2B  coSfficienten  bepaald  uitgaande  van  verschillende  startwaarden 


n 

h 

s* 

•o 

A 

opt-'i 

Opt_c^ 

opt_®o 

Opt_A 

«c  lt«r 

1 

0.00019 

0.00010 

0.01900 

46.063 

24.1630 

1 

0.00010 

0.00010 

0.00001 

10.000 

13.4146 

9 

2 

0.00019 

0.00010 

0.19000 

46.063 

16.0932 

1 

0.00010 

0.00000 

0.00002 

10.000 

13.4022 

3 

100* 

3 

0.00019 

0.00010 

1.90000 

46.063 

4.3311 

1 

0.00192 

0.00115 

0.19036 

46.494 

0.3795 

0 

19 

4 

0.00019 

0.00100 

0.01900 

46.063 

5.4740 

1 

0.00192 

0.00115 

0.19036 

46.494 

0.3795 

0 

24 

S 

0.00019 

0.00100 

0.19000 

46.063 

3.0749 

1 

0.00192 

0.00115 

0.19036 

46.494 

0.3795 

0 

19 

6 

0.00019 

0.00100 

1.90000 

46.063 

7.0366 

1 

0.00192 

0.00115 

0.19036 

46.494 

0.3795 

0 

25 

7 

0.00019 

0.01000 

0.01900 

46.063 

15.4402 

1 

0.00192 

0.00115 

0.19036 

46.494 

0.3795 

9 

32 

1 

0.00019 

0.01000 

0.19000 

46.063 

15.4966 

1 

0.00192 

0.00115 

0.19036 

46.494 

0.3795 

0 

34 

9 

0.00019 

0.01000 

1.90000 

46.063 

15.9790 

1 

0.00192 

0.00115 

0.19036 

46.494 

0.3795 

0 

31 

10 

0.00190 

0.00010 

0.01900 

46.063 

24.1630 

1 

0.00190 

0.00003 

0.01653 

10.000 

13.4090 

9 

n 

0.00190 

0.00010 

0.19000 

46.063 

IS. 7500 

1 

0.00010 

0.00003 

0.00026 

35.746 

3.7313 

3 

100* 

12 

0.00190 

0.00010 

1.90000 

46.063 

4.6596 

1 

0.00192 

0.00115 

0.19036 

46.494 

0.3795 

0 

14 

13 

0.00190 

0.00100 

0.01900 

46.063 

3.7156 

1 

0.00192 

0.00115 

0.19036 

46.494 

0.3795 

0 

29 

14 

0.00190 

0.00100 

0.19000 

46.063 

1.4366 

1 

0.00192 

0.00115 

0.19036 

46.494 

0.3795 

0 

16 

15 

0.00190 

0.00100 

1.90000 

46.063 

7.6369 

1 

0.00192 

0.00115 

0.19036 

46.494 

0.3795 

0 

23 

1< 

0.00190 

0.01000 

0.01900 

46.063 

15.6934 

1 

0.00192 

0.00115 

0.19036 

46.494 

0.3795 

0 

40 

17 

0.00190 

0.01000 

0.19000 

46.063 

15.7506 

1 

0.00192 

0.00115 

0.19036 

46.494 

0.3795 

0 

34 

11 

0.00190 

0.01000 

1.90000 

46.063 

26.2214 

1 

0.00192 

0. 00115 

0.19036 

46.494 

0.3795 

0 

31 

19 

0.01900 

0.00010 

0.01900 

46.063 

9.9651 

1 

0.00192 

0.00115 

0.19036 

46.494 

0.3795 

0 

29 

20 

0.01900 

0.00010 

0.19000 

46.063 

12.7464 

1 

0.00192 

0.00115 

0.19036 

46.494 

0.3795 

0 

25 

21 

0.01900 

0.00010 

1.90000 

46.063 

304.0929 

1 

0,00192 

0.00115 

0.19036 

46.494 

0.3795 

0 

49 

22 

0.01900 

0.00100 

0.01900 

46.063 

15.7194 

1 

0.00010 

0.00150 

0.00001 

45.979 

0.5407 

9 

_ 9 

23 

0.01900 

0.00100 

0.19000 

46.063 

17.1024 

1 

0.00192 

0.00115 

0.19036 

46.494 

0.3795 

0 

32 

24 

0.01900 

0.00100 

1.90000 

46.063 

249.6360 

{ 

0.04744 

0.00000 

7.01921 

10.000 

99.1024 

0 

22* 

25 

0.01900 

0.01000 

0.01900 

46.063 

16.3241 

» 

0.00192 

0.00115 

0.19036 

46.494 

0.3795 

0 

41 

26 

0.01900 

0.01000 

0.19000 

46.063 

16.3756 

1 

0.00192 

0.00115 

0.19036 

46.494 

0.3795 

0 

39 

27 

0.01900 

0.01000 

1.90000 

46.063 

16.6174 

1 

0.00192 

0.00115 

0.19036 

46.494 

0.3795 

0 

46 

Commentaar 

De  grafieken  tonen  dat  de  origineel-  en  resultaatpolylijn  niet  perfect  op  elkaar  liggen. 
Voor  de  kleinste  Opt.A  bedraagt  de  veitiouding  Opt_A/Opt_0k)  ongeveer.  0.8  %. 
(0.3795  /  46.494) 

AUe  dechie  tupels  onderscheiden  zich  duidelijk  door  de  sterk  afwi^cende  OpCA. 

De  enige  matige  ttq>el  (22)  kan  op  fysisdK  gnmden  uitgesloten  woiden  aangezien  Opt_c« 
een  randwaaide  van  aameemt  van  het  zodtgebied. 

Gevcdg:  er  reafeien  enkel  foede  tupds  en  er  is  6dn  uniek  mintmum  gevondea 

Het  schattenlgoritine  (met  ststwaaidetdjqMling)  heeft  het  atwdide  minimmn  in  het 

zoekfdjied  fevonden. 
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Tabel  5.  Land  coeflicienten  uitgaande  van  versdiillende  startwaarden 


n 

'2 

Si 

•o 

A 

Opt_f2 

Opt_C^ 

Opt_«lo 

Opt^A 

VC  It«r 

1 

0.00130 

0.00010 

0.01900 

€3.471 

3€.442€  1 

0.00130 

0.00004 

0.00009 

55. €24 

2.ffl€ 

3 

100* 

2 

0.00130 

0.00010 

0.19000 

€3.471 

34.4142  1 

0.00010 

0.00001 

0.00001 

55.500 

2.920f 

3 

100* 

3 

0.00130 

0.00010 

1.90000 

€3.471 

€.1304  1 

0.01349 

0.00101 

0.19409 

€3.440 

0.1050 

0 

23 

4 

0.00130 

0.00100 

0.01900 

€3.471 

17.243S  1 

0.01349 

0.00101 

0.19409 

€3.440 

0.1050 

0 

23 

S 

0.00130 

O.OQIOQ 

0.19000 

€3.471 

11.4017  1 

0.01349 

0.00101 

0.19409 

€3.440 

0.1050 

0 

20 

6 

0.00130 

0.00100 

1.90000 

€3.471 

10.237S  1 

0.01349 

0.00101 

0.19409 

€3.440 

0.1050 

0 

34 

7 

0.00130 

0.01000 

0.01900 

€3.471 

20.3552  1 

0.01349 

0.00101 

0.19409 

€3.440 

0.1050 

0 

44 

1 

0.00130 

0.01000 

0.19000 

€3.471 

20.4175  1 

0.01349 

0.00101 

0.19409 

€3.440 

0.1050 

0 

4€ 

9 

0.00130 

0.01000 

1.90000 

€3.471 

22.51f€  1 

0.01349 

0.00101 

0.19409 

€3.440 

0.1050 

0 

4f 

10 

0.01300 

0.00010 

0.01900 

€3.471 

20.0253  1 

0.01349 

0.00101 

0.19409 

€3.440 

0.1050 

0 

23 

u 

0.01300 

0.00010 

0.19000 

€3.471 

11.7972  1 

0.01349 

0.00101 

0.19409 

€3.440 

0.1050 

0 

24 

12 

0.01300 

0.00010 

1.90000 

€3.471 

12.7259  ( 

0.01349 

0.00101 

0.19409 

€3.440 

0.1050 

0 

20 

13 

0.01300 

0.00100 

0.01900 

€3.471 

4.9925  1 

0.01349 

0.00101 

0.19409 

€3.440 

0.1050 

0 

17 

14 

0.01300 

0.00100 

0.19000 

€3.471 

0.4392  1 

0.01349 

0.00101 

0.19409 

€3.440 

0.1050 

0 

7 

13 

0.01300 

0.00100 

1.90000 

€3.471 

15.4722  1 

0.01349 

0.00101 

0.19409 

€3.440 

0.1050 

0 

40 

1( 

0.013C0 

0.01000 

0.01900 

€3.471 

23.1023  1 

0.01349 

0.00101 

0.19409 

€3.440 

0.1050 

0 

3€ 

17 

0.01300 

0.01000 

0.19000 

€3.471 

23.3332  1 

0.01349 

0.00101 

0.19409 

€3.440 

0.1050 

0 

39 

If 

0.01300 

0.01000 

1.90000 

€3.471 

35.0441  ) 

0.01349 

0.00101 

0.19409 

€3.440 

0.1050 

0 

45 

19 

0.13000 

0.00010 

0.01900 

€3.471 

149.0012  1 

0.01349 

0.00101 

0.19409 

€3.440 

0.1050 

0 

37 

20 

0.13000 

0.00010 

0.19000 

€3.471 

249.4429  1 

0.01349 

0.00101 

0.19409 

€3.440 

0.1050 

0 

33 

21 

0.13000 

0.00010 

1.90000 

€3.471 

53S.3C94  1 

0.01349 

0.00101 

0.19409 

€3.440 

0.1050 

0 

59 

22 

0.13000 

0.00100 

0.01900 

€3.471 

125.0520  i 

0.01349 

0.00101 

0.19409 

€3.440 

0.1050 

0 

39 

23 

0.13000 

0.00100 

0.19000 

€3.471 

179.1244  1 

0.01349 

0.00101 

0.19409 

€3.440 

0.1050 

0 

31 

24 

0.13000 

0.00100 

1.90000 

€3.471 

913.0240  1 

0.01349 

0.00101 

0.19409 

€3.440 

0.1050 

0 

34 

25 

0.13000 

0.01000 

0.01900 

€3.471 

53.5453  1 

0.01349 

0.00101 

0.19409 

€3.440 

0.1050 

0 

52 

2« 

0.13000 

0.01000 

0.19000 

€3.471 

54.5294  1 

0.01349 

0.00101 

0.19409 

€3.440 

0.1050 

0 

€0 

27 

0.13000 

0.01000 

1.90000 

€3.471 

914. 9444  1 

0.00010 

0.02177 

0.00001 

100.000 

24.2750 

3 

100* 

Optnerkingen: 

De  grafieken  tonen  dat  de  origineel-  en  resultaaqpolylijn  keurig  op  elkaar  liggen. 

Voor  de  kleinste  Opt_A  bedraagt  de  vertiouding  Opt_A/C)pt_(ii)D  ongeveer  0.2  %. 
(0.1050  /  63.440) 

Alle  slechte  tupels  onderscheiden  zidi  duidelijk  door  de  steik  afwijkende  OpLA  Gevolg: 
er  resteten  enkel  goede  tupels  en  er  is  ddn  uniek  minimum  gevondea 
Het  schatteralgoritme  (met  startwaardenbepaling)  heeft  het  abadtde  minimum  in  het 
zoekgebied  gevondea 
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10^  Testconclusies 


De  grafieken  van  UitrollA,  UitiDl2A  en  Land  tonen  dat  de  origineel-  en  de  resultaaqMlylijnoi 
keurig  op  elkaar  liggen.  Bij  UitrollB  en  Uitrol2B  is  dit  niet  het  geval.  De  coefBcienten  van 
UitrollB  en  Uitrol2B  zijn  deituilve  niet  betrouwbaar. 

Voor  de  kleinste  Opt_A  ligt  de  vertiouding  Opt_A  /  Opt_aio  steeds  tussen  0.2  %  en  1.3  %,  waaibij 
de  slechtste  (hoogste)  veihoudingen  afkomstig  zijn  van  de  B-polylijnen.  Voor  een  optisch  'goede* 
uitrolcurve  bedraagt  de  vertiouding  Opt_A  /  Opt_(0o  hoogstens  0.5  %. 

Ofschoon  niet  gegarandeerd  kan  worden  dat  het  schatteralgoritme  het  beste  minimum  binnen  het 
zoekgebied  vindt,  heeft  het  algoritme  in  alle  geteste  gevallen  dit  beste  minimum  blijkbaar  wel 
gevonden. 

Alle  slechte  tupels  onderscheiden  zich  duidelijk  door  de  stertc  afWijkende  Opt_A. 

Bij  matige  tupels  dient  gekeken  te  wonten  naar  de  waarden  van  de  fysische  co^fTicienten.  Indien 
ddn  of  meer  coiffficiSnten  een  randwaarde  van  het  zoekgebied  aannemen,  dan  zijn  de  cotifficienten 
niet  betrouwbaar. 

Het  is  overigens  onwaarschijnlijk  dat  het  schatteralgoritme  leidt  tot  een  matige  tupel  zoals  bij 
UitrollB,  Uitrol2A  en  UiUol2B  is  vooigekomen.  De  staitwaarden  zoals  berekend  m.b.v.  de 
methode  der  kleinste-kwadraten  blijken  namelijk  reeds  erg  goed  te  zijn: 


Tabel  6.  Schatteralgoritme 


b«ttand 

'2 

Sr 

•o 

A 

UitrollA 

0.00<70 

0.00040 

0.37111 

40.270 

0.2007 

UitrollB 

0.00700 

0.00020 

0. 40551 

40.270 

0.7124 

Uitrol2A 

0.00750 

0.00033 

0.42000 

40.003 

0.1755 

ultrolBB 

0.00101 

0.00134 

0.04551 

40.003 

0.4200 

LoftB 

0.01235 

0.00124 

0.13700 

03.471 

0.1005 

apt-«2 

opt.c. 

Opt.wo 

Opt_A 

«C  xt«r 

0.00700 

0.00031 

0.43555 

40.350 

0.2000 

0  0 

0.00750 

0.00000 

0.70042 

40.702 

0.5044 

0  0 

0.00773 

0.00031 

0.42722 

45.030 

0.1030 

0  7 

0.00102 

0.00115 

0.10030 

40.404 

0.3705 

0  11 

0.01340 

0.00101 

0.10400 

03.440 

0.1050 

0  0 
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De  A's  behorende  bij  de  matige  tupels  zijn  ongeveer  een  factor  100  groter  dan  de  A's  afkomstig 
van  de  met  bet  schatteralgoritme  berekende  startwaarden.  De  A's  behorende  bij  de  staitwaarden 
zoals  berekend  door  het  schatteralgoritme  liggen  al  redelijk  dicht  bij  Opt_A.  Dit  geeft  aan  dat  de 
kleinste-kwadraten  startwaaidenbepaling  goed  voldoet. 

Omdat  zowel  UitiollA  als  Uitrol2A  goede  uitrolcurves  zijn  en  beide  onder  dezelfde 
omstandigheden  en  met  hetzelfde  voertuig  zijn  uitgevoeid,  is  het  zinvol  de  fysische  co£ffici£nten 
met  elkaar  te  vergelijken  (zie  tabel  6).  De  fysische  cotffficiSnten  behoren  ongeacht  de  uitiolcurve 
identiek  te  zijn.  Relatief  beschouwd: 

Opt_fi(UitiollA)  -  Opt_fi(Uitrol2A) 

-  =-1.7% 

Opt_fi(UitrollA) 

Opt_f2(UitrollA)  -  Opt_f2(Uitrol2A) 

-  =  0.0  % 

Opt_f2{UitrollA) 

Opt_c^(Uitroll  A)  -  Opt_c^(Uitrol2A) 

-  =1.9% 

Opt_c»(UitrollA) 

Op  grond  van  deze  veitioudingen  kan  niet  beweerd  woiden  dat  de  coSflicienten  met  een  precisie 
van  2  %  bepaald  kuimen  woiden.  Hiervoor  zijn  er  te  weinig  A-uitrolcurves. 

Merk  op  dat  het  niet  zinvol  is  een  soortgelijke  verhouding  voor  Opt.capCUitiollA)  en 
Op(_oiQGJitit)I2A)  te  berekenen.  (Oq  staat  immers  voor  de  snelheid  ten  tijde  van  de  eerste  meting 
en  is  dus  altijd  proefafhankelijk. 

Het  ingestelde  maximum  aantal  iteraties  blijkt  voldoende  groot  te  zijn.  opdat  het  schatteralgoritme 
een  minimum  ktui  vinden. 

T.a.v.  de  waarde  van  dient  opgemerkt  te  woiden  dm  de  foimule  ingebouwd  in  het 
schatteralgoritme  waar  de  testen  mee  uitgevoeid  zijn  een  afleveit  die  factor  2  te  klein  is.  De 
c^-waaide  van  de  DAF  wordt  na  coirecde  ongeveer  0.86.  Deze  waarde  ligt  in  het  interval  van 
[0.8;  2.0]  hetgeen  door  KraftfahH<^  opgegeven  wordt  als  een  normde  waarde  voor  vradttwagens. 
De  c^-waarde  van  0.38  voor  de  landrover  is  niet  geverifieerd  met  fabridcsgegevens. 


De  veihouding  Opt_fi  /  Opt_f2  is  met  een  ordc  van  gnx>tte  10  veel  kleiner  dan  op  dc  theoretische 
gronden  verwachte  waarde  van  de  oide  van  gnootte  1000.  Hiervoor  is  geen  veiUaring  gevonden. 

De  grillige  A-versdiilcurves  wekken  niet  de  indruk  dat  er  belangrijke  effecten  ajn  die  door  het 
fysische  model  niet  beschreven  woidea  Het  meenemen  van  de  windsnelheid  in  het  fysisch  model 
is  echter  het  overwegen  waard.  Dit  zal  vermoedelijk  slechts  een  gering  effect  hebben  op  Opt_A, 
maar  heeft  mogelijk  anderc  waaiden  voor  de  fysische  codfficienten  tot  gevolg. 

Merk  op  dat  Warning  Code  8  (=  Mogelijk  niet  optimale  minimum  gevonden),  afkomstig  van 
DBCONF,  niet  gebruikt  is  om  tupels  af  te  wijzen. 


10.6  GEBRUIKERSADVIEZEN 


De  resultaten  van  de  testen  van  het  schatteraigoritme  leiden  tot  de  volgende  gebniikersadviezen: 
Het  veidient  de  aanbeveling  meer  dan  66n  uitrolcurve  te  maken  met  hetzelfde  voettuig 
onder  dezelfde  omstandigheden. 

Er  diem  gecontroleerd  te  worden  of  de  co^fficiCnten  randwaaiden  van  het  zoekgebied 
aannemen.  Indien  dit  het  geval  is,  zijn  de  co€ffici(inten  waarschijnlijk  onbetrouwbaar. 

Er  dient  zo  mogelijk  gebruik  gemaakt  te  worden  van  gefixeeide  codfficientea 
Er  dienen  geen  veikorte  uitrolcurves  gebruikt  te  worden.  Sc  bijlage  C. 
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11  CONCLUSIE 


Het  schatteralgoritme  zoals  gepresenteerd  in  dit  document  is  in  staat  wnjvingscot5ffici6nten  te 
berekenen  voor  de  UDB.  Om  tot  statistische  (algcmcne)  uitspraken  te  kcanen  m.b.t  tot  de 
nauwkeurigheid  van  het  schatteralgoritme  is  een  veel  groter  aantal  goede  uitrolcurves 
noodzakelijk.  De  rcsultaten  zoals  door  het  schatteralgoritme  geleverd  voor  de  uitrolcurves 
UitrollA  en  UitroUA  zijn  mathematisch  bcschouwd  heel  bevrcdigend  vanwcge  de  in  absolute  zin 
lage  nonn.  De  overeenstemming  is  fysisch  gezien  heel  hoopgevend.  De  berekende  c„*waarde  van 
de  DAF  YA-4440  komt  overeen  met  hetgeen  volgens  de  literatuur  verwacht  mag  worden.  De 
wrijvingsco£ffici£nten  Opt_fj  en  Opt_f2  voldoen  niet  helemaal  aan  de  fysische  verwachtingen.  De 
reden  hiervoor  is  onbekend.  Of  de  coCfficignten  in  de  prakUjk  bevrcdigend  zijn.  kan  pas  na 
gebruik  in  de  rijsimulatorgegvalueerd  worden. 
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DE  METHODE  VAN  EULER 


Een  stelsel  differentiaalvergelijkingen  kan  m.b.v.  de  methode  van  numeriek  benadeid 

worden. 

Zij  gegeven  de  volgende  differentiaalvergelijking  m6t  randvoorwaarde: 
jU’(t)  =f(t.U(t)) 

^  U(to)  =  Co 

waaxbij  f  e  C[cx,p]  *  D  (=  Verzameling  continue  fiincties  op  [a,p]  *  D.  en  D  c  R),  Cq  e  R.  Kies 
stapgrootte  h  >  0  en  definieer:  t,  =  a  +  n  *  h  (n  =  0,  1,  2, ...).  Zij  N  het  gnx)tste  gehele  getal  z6 
dat  In  S  p.  Dan  is  Un  een  benadering  voor  U(t^  met: 
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DE  METHODE  DER  KLEINSTE  KWADRATEN 

Zij  gegeven  {(Xj,  yj)  e  R2  I  i  =  1,  2 . N).  Gevraagd  dc  codfficiaiten  U],  02  en  (X3  van  het 

polynoom: 

y(x)  =  a,  +  02  X  +  03  x2 
zd  dat 

N 

I  (y(Xi)-yi)2 

i=l 

minimaal  is. 


Zij  A  3:  {ajj}  een  N*3  matrix,  z6  dat: 


ai,  =  l 

i=1.2 . N 

ai2  =  Xi 

i=  1,2,  ...,N 

aa  =  Xi2 

i=  1,2,  ...,N 

Laat  A'^  de  getransformeerde  van  A  zijn.  Dan  worden  de  klcinste  kwadraten  codfTiciCnten  O],  02 
en  03  bepaald  door  het  volgende  stelsel  lineaire  vergelijkingen  op  te  lossen: 

(oil)  (yi) 

ATA(a2)  =AT(...) 

(013)  (yN) 


HEX  GEBRUnC  VAN  VERKORTE  UTTROLCURVES 


In  deze  bijlage  woidt  ingegaan  op  de  vraag  of  het  mogelijk  is  de  snelheid  als  fiinctie  van  de  tijd 
met  een  gedeelte  van  de  uitrolcurve  te  besdirijven.  Indien  een  voertuig  niet  behoeft  te  accelereien 
tot  een  hoge  aanvangssnelheid  of  uit  te  rollen  tot  stilstand.  kan  namelijk  volstaan  worden  met  een 
koitere  testweg.  Dit  komt  de  veikeersveiligheid  ten  goede  en  vergroot  het  aantal  geschikte 
testwegen.  Geschikte  testwegen  voor  het  maken  van  uitrolcurves  van  zware  voertuigen  zijn 
schaars  aangezien  deze  zo'n  2  km  lang  en  recht  dienen  zijn. 

Helaas  blijkt  het  niet  mogelijk  te  zijn  om  de  uitrolcurves  te  verkorten.  Daartoe  een  vooibeeld. 

In  figuur  6  zijn  twee  polylijnen  goed  te  onderscheiden.  Als  veikoite  uitrolcurve  is  uitrolcurve  lA 
genomen.  met  uitsluitend  die  metingen  waarvan  de  gemeten  snelheden  in  het  interval 
[30  km  /  h,  70  km  /  h]  liggen.  De  naar  boven  afbuigende  lesultaa^iolylijn  in  figuur  6  is  berekend 
met  de  coefficien  :n  die  optimaal  zijn  t.o.v.  van  de  veikone  uitrolcurve.  De  andeie  polylijn  is  de 
origineelpolylijn  van  uitrolcurve  1 A  met  een  snelheidsbereik  van  0  tot  70  km  /  h. 
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Fig.  6.  Verkorte  lA-uitrolcurve. 


Uit  de  grafiek  volgt  dat  de  coSHiciSnten  berekend  met  een  veikorte  uitrolcurve  niet  de  gehele 
uitrolcurve  beschrijven,  tioewel  ze  wel  een  goede  fit  leveren  op  het  beperkte  stuk. 
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